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Preface 



The star formation (SF) can be regarded as an engine powering the light emission 
of spiral galaxies at several wavelengths. Massive, young stars emit the bulk of 
their light in the UV band; a large fraction of the UV radiation is absorbed by cold 
(T ~ 40K) dust and reradiated in the far infrared band. Massive stars explode as 
supernovae, and supernova remnants accelerate the electrons which in turn power 
the radio emission through the synchrotron mechanism. Binary systems composed 
by a collapsed star with a massive companion, emit X-rays through Comptonization 
of thermal photons coming from the accretion disc. 

The star formation activity may go on for very long times at low paces, or be 
significantly enhanced for short times as conquence of perturbations in the gravity 
potential, which trigger the collapse of gas clouds. Mergers, or cannibalization of 
small galaxies, are often the reason for episodes of enhanced star formation. 

The intensity of an episode of star formation can be classified on the basis of its 
sustainability for long times. If we define the Star Formation Rate (SFR, usually 
expressed in M /yr) as the mass of gas converted into stars per unit time, we may 
also consider the depletion time r d = M gas /SFR which measures the time needed to 
convert into stars all the mass M gas of gas present in a galaxy if the star formation 
continues at a constant rate. If the depletion time is very short (much shorter than, 
say, a Hubble time) we may call the present SF episode a starburst episode. The 
galaxies with an undergoing starburst episode can be called starburst galaxies. Note, 
however, that in the literature the term "starburst galaxy" is also found to be loosely 
referred to large galaxies with quiescent SF but large SFR values. 

The star formation activity plays a significant role in the evolution of galaxies. 
The production of heavy chemical elements results both in the build up of different 
stellar populations and in the pollution of the insterstellar medium (ISM). Intense 
starburst episodes also may form bubbles of hot ionized gas, with a pressure larger 
than that of the surrounding cold ISM. The bubbles may be able to expand in a 
direction normal to the plane of the galaxy disc, eventually releasing their gas content 
into the intergalactic medium (IGM). Depending on the strength of gravitational 
potential of the galaxy, the gas may definitely escape into the IGM or fall back (a 
galactic fountain is then formed). 

Star formation has also undergone a strong cosmological evolution, since the av- 
erage density of SFR was about ten times larger at redshifts la rger than 1, and has 
been declining since then ( Lilly et al. 19961: Madau et al. 19961 ). Thus the study of 
galaxies in the local universe with a strong star formation activity may act as a 
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Preface 



guide to understand the evolution of galaxies in the early universe. 

In the following, we focus on the study of 'normal', spiral, star forming galaxies. 
The term 'normal' will be referred to galaxies whose energy output is dominated 
by all the processes related with star formation and evolution, with no contribution 
from a possible active galactic nucleus. 

The first chapter is devoted to a review of the infrared, radio and X-ray emission 
from star forming galaxies. In the second chapter we explore the quantitative re- 
lationships between the total X-ray, radio and infrared luminosities for a sample of 
23 star forming galaxies. It is found that linear correlations hold between the lumi- 
nosities in the above bands. In the third chapter we show, by analyzing a sample 
of 11 high redshift galaxies in the Hubble Deep Field, that the linear correlations 
may be extended up to z ~ 1. In the fourth chapter we turn to an analysis of the 
luminosity function and number counts of star forming galaxies; we show that the 
number density of normal galaxies at faint X-ray fluxes (10 _17 -10 -15 erg s _1 cm~ 2 ) 
is very well defined. 

The fifth chapter offers a first insight into a different aspect of galaxy evolution, 
i.e. whether is it possible to determine the metallicity enhancements in a starburst 
episode. An interesting, yet problematic, picture of metal enrichment of the inter- 
stellar matter arises from an analysis of the X-ray and near infrared spectra of M82. 
Finally, a summary is offered of the work here described. 
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1 Light emission from star forming 
galaxies 



Two main stellar population are usually found in spiral galaxies: one made up by 
old dwarf stars, and one made up by young, blue, massive stars. While the former 
builds the bulk of the galaxy stellar mass and of the energetic output in the optical 
band, the latter powers the light emission in many different bands throughout the 
whole energy spectrum. In the following, it will be shown how the processes of 
star formation and death can be invoked to explain the light emission properties of 
galaxies at radio, far infrared and X-ray wavelengths. 



1.1 Infrared emission from spiral galaxies 

The far infrared luminosity (~ 3-300/i) of a star forming galaxy is commonly inter- 
preted as due to starlight reradiation by dust grains. The early evid ence for a link be- 



tween infrared emission and dust was based on several arguments (|Harwit fc Pacini 
197,4 iTelesco fc Harper! Il980l : iDevereux fc Younglll990h : 



• H II regions are associated with knots of infrared emission; 

• the far infrared luminosities may exceed the optical luminosities, consistent 
with the view that the FIR radiation is generated by heavily obscured star 
clusters; 

• the high luminosity-to-mass ratios derived from the comparison of FIR lumi- 
nosities and dynamical masses, consistent with the view that the obscured star 
clusers are composed primarily of young massive stars; 

• the FIR luminosities are well correlated with masses of the molecular gas, 
which is the component of the interstellar medium from which stars form. 

The observations made with the Infrared Astronomical Satellite (IRAS) showed that 
infrared emission is ubiquitous among spiral galaxies, thus reinforcing the view that 
the far infrared luminosity of spiral galaxies is thermal dust emis sion. A comparison 



of th e Ha and FIR luminosities in a sample of 124 galaxies ([Devereux fc Young 



1990) showed that the mean ratio of extinction-corrected Ha to bolometric FIR 



luminosity in spiral galaxi es is comparable to that ex pected for H II regions powered 
by massive stars. This led Devereux fc Youriel (|l990h to support the view that high 
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1 Light emission from star forming galaxies 
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Figure 1.1: Far infrared spectral en- 
ergy distributions for a sample of spi- 
ral ga laxies in Virgo, from Popescu et al.l 
(|2002h . The points are ISOPHOT data, 
while the lines show the two-temperature 
model fit to the data (see text). 
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Figure 1.2: Fractional contribution 
of ultraviolet light from different sources 
(Hn regions and the diffuse old stellar 
component) to dust hea ting for the spiral 
galaxy NGC 891, from iTuffs fc Popescu 



mass stars required to ionize the hydrogen gas can easily generate the far infrared 
luminosities. 

However, the presence of galaxies with a low level of i nfrared emission that is ener- 
getically comparable to the optical luminosity also led lde Jong et al. (|l984h to sug- 
gest later type stars as a significant dust heating source, p articularly in spirals with a 
cool fffinw /ffinn/, dust tempera ture. It was also suggested (jLonsdale-Persson fc Helou 
19871 : Ide Jong fc Brinklll987l ) that the emission at and 100/i measured by IRAS 
could be deconvolved in two components: a warm component with temperatures 
of 50-60K, and a cool component (the cirrus) at 10-20K. An interpretation of the 
two-component model which considered the warm component as radiation from dust 
heated by young stars, while the cool component is heated by older, late type stars, 
was conside red as the possib le explanation for a slight non-linearity of the radio/FIR 
correlati on dFitt et all 19881 however, the latest studies point toward almost exact 
linearity, Yun et al. l200ll : see Sect. 11.3ft . 

Since the peak in the thermal spectrum of dust with temperatures below ~ 30K 
falls longwards of 100/i — thus outside the IRAS bands — , a clear separation of the 
warm and cool components was only possible with the launch of the Infrared Space 
Observatory (ISO), w hose spectral coverage extended up to wavelengths of ~ 180/x. 
Popescu et al. (2002) studied a sample of 38 galaxies with the ISOPHOT instrument 
at 60, 100 and 170/a. Most of these galaxies were discovered to contain a cold dust 
emission component which could not have been detected by IRAS. The FIR Spectral 
Energy Distributions (SEDs) of these objects were fitted with a combination of two 
modified blackbody functions, one physically identified with the emissions from a 
localised warm dust component, with a fixed temperature of 47K, and one other 
with a diffuse cold dust component (Figs. 11.11 11.2ft . The cold dust temperatures 
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1.1 Infrared emission from spiral galaxies 




Fig ure 1.3: ISOPHOT mapping of the M31 galaxy at 170/i, from lHaas et al. 
The darker areas show the emission due to the cold (~ 20 K) dust. The 
angular resolution is 1.3', the field size 2.9° x 2.9°. 
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1 Light emission from star forming galaxies 

were found to span a broad range, with a median of 18K. 

Imaging of the cold component has been so far possible, with ISOPHOT, only for 
very extended objects such as M31 (Fig. I1..3J) . Detailed studies of the FIR emission 
will be routinely possible in the immediate future with the Spitzer Space Telescope. 
Among the first images taken with this new observatory, the pictures of M81 at 
different wavelengths ( Plates HUT]) show the details in the spiral structure of the cold 
dust at 170//. 



1.2 Radio emission 

Radio sources exist in nearly all normal spiral and dwarf irregular galaxies, and in 
many peculiar or interacti ng systems. Th e morphology of the radio emission found in 
normal galaxies is various (jCondonl 1992, and references therein) . The faintest dwarf 
irregulars have radio luminosities comparable with the Galactic Supernova Remnant 
(SNR) Cassiopea A and little or no detectable emission extending beyond known H II 
regions and SNRs. Their radio morphologies are lumpy and irregular (see IC 10 in 
Fig. ll.4|) . and their radio spectra are often relatively flat. The thick radio disk/halo of 
the edge-on galaxy NGC891 and the fairly smooth radio disk with bright spiral arms 
of the face-on galaxy NGC6946 are typical of the larger spiral galaxies. The central 
radio sources in normal galaxies like NGC6946 have a larger brightness but are 
usually much less luminous than the disk-like sources. Luminous radio sources with 
complex morphologies may be found in colliding galaxies (e.g. NGC1144). There is 
a tendency for fairly compact (diameter D < 1 kpc) central starbursts to dominate 
at higher radio luminosities, as in M82. The most luminous radio sources in normal 
galaxies are frequently quite compact (D ~ 200 pc), and confined to the nuclei of 
strongly interacting systems (e.g. IC694+NGC3690). Numerical simulations show 
that a collision involving a disk galax y can drive about half of the disk gas (< 



1O 1O M ) within 200 pc of the nucleus (jHernquistl 119891 : iBarnes Hernquistl ll991 ) 



where the gas de nsity becomes quite high before a powerful starburst is triggered 
( KennicuttllSl . The resulting massive stars and their SNRs then produce intense 



radio emission. 

Nearly all of the radio emission from normal galaxies is synchr otron radiatio n 
from relativistic electrons and free-free emission from Hn regions (|Gondonlll992h . 



Thermal reradiation of starlight by dust quickly overwhelms these components above 
v ~ 200GHz, defining a practical separation between the radio and infrared bands. 
Typical intensities of synchrotron radiation, free-free emission, and dust reradiation 
are shown in the radio/FIR spectrum of M82 (Fig. II the radio continuum only 
accounts for < 10~ 4 of its bolometric luminosity. 

However, determining the amount of the thermal fraction is highly controversial. 
In principle, the relatively flat spectrum (S(u) oc v~ a with a < 0.1) thermal emission 
should be distinguishable from the steeper spectrum [a ~ 0.8) nonthermal emission 
via total flux densities or maps obtained at two or more frequencies. In practice, 
most normal galaxies are not bright enough to be detected at frequencies much 
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1.2 Radio emission 









Figure 1.4: Contour maps illustrating the range of radio source morphologies, 
sizes, and luminosities found in normal galaxies, from Condon! (|l992| ). The bars are 
2/i -1 kpc long. The logarithmic contours are separated by v^2 in brightness, and 
the 1.49 GHz brightness temperatures of the lowest contours are 0.25K (IC10, 
NGC891, NGC6946), 0.5K (M82), 8K(NGC1144), and 128K (IC694+NGC3690). 



higher than v ~ 10 GHz (A ~ 3 cm). Moreover, if integrated flux densities at only 
three or four frequencies are used to fit the (unknown) nonthermal spectral index and 
the thermal emission simultaneously, the resulting nonthermal spectral index and 
thermal fraction are strongly correl ated. An upper limit St I S < 0.4 to the thermal- 



to-total flux ratio was estimated by lKlein fc Emerson! ()198l! ) and lGioia et al.l (Jl 



from the lack of a flattening of the radio spectrum in two samples of nearby spiral 
galaxies. The thermal fraction is highly uncertain even in the brightest and best 
observed galaxies such as M82 because the nonthermal compon ent may not have a 



straight spectrum. By assuming a power law spectrum for M82. iKlein et al. 
calculated a thermal flux density St = 0.15 Jy at v = 32 GHz. Using practically the 
same data but assuming that the nonthermal spectrum steepens at high frequencies, 
Carlstrom fc Kronberd (1991) obtained St = 0.5 Jy (and thus a thermal fraction of 
~ 70%, cfr. Fig. ESI) at v = 92 GHz, where S T = 0.17 Jy (thermal fraction ~ 25%) 
would be expected by extrapolating Klein's estimate at 92GHz with a typical slope 
a = 0.1. This simple consideration shows that the estimates of the thermal fraction 
at high frequencies may have un uncertainty of (at least, since M82 is among the 
best studied galaxies) a factor of three. 
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1 Light emission from star forming galaxies 
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Figure 1.5: Radio emission from M82. Upper panel: medium r esolution image 
at 5GHz (Cambridge 5km radio telescope) . from lHargrave (Il974 ) . Lower panel: 
high resolution Merlin image at 5GHz, from iMuxlow et aL ( 1994f ). As normal with 
interferometer images, the lower the resolution, the higher the sensitivity to diffuse 



emission. 
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1.2 Radio emission 
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Figure 1.6: Radio-FIR spectrum of M82, from iCondonl (jl992f ). Dashed line: 
thermal contribution to the spectrum from the Hll regions. Dot-dashed line: non- 
thermal emission. Dotted line: black body emission from the dust. The continous 
line shows the fitted spectrum resulting from the sum of the previous three compo- 
nents. 



Since the nonthermal spectrum has a steeper slope, the uncertainty on the ther- 
mal fraction becomes a minor issue at lower frequencies (e.g. at 1.4 GHz), where the 
majority of the radio emission is of nonthermal origin (cfr. Fig. ll.fiJ) . The nonthermal 
radiation is usually explained as synchrotron radiation from relativistic electrons. 
Only stars more massive than ~ 8M produce the core-collapse supernovae whose 
remnants (SNRs) are thought to accelerate most of the relativistic electrons in nor- 
mal galaxies; these stars ionize the Hll regions as well. Their lifetime is < 3 ■ 10 7 yr, 
while the relativistic electrons probably have lifetimes < 10 8 yr. Radio observations 
are therefore probes of very recent star formation activity in normal galaxies. 

Supernova remnants become radio sources about 50 years after the explosion, as 
Rayleigh- Taylor instabilities develop in the boundary between the shock and the am- 
bient interstellar medium, and remain visibile for hundreds or thousands of years. 
About 40 young SNRs are conspicuos in high resolution maps of M82 at 5 GHz 
f)Kronberg et al.lll985l : iMuxlow et al.lll994L Fig. EH). Although SNRs are probably 
responsible for cosmic ray acceleration, discret e remnants themselves emit only a 
small fraction (< 10%) of the integrate d flux (|Poolev 19691 : Ilovaiskv fc Lequeux 
19721: lBiermannlll976l: i Helou et al.lll985h . This also means that > 90% of the non- 



thermal emission must be produced long after the individual supernova remnants 
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1 Light emission from star forming galaxies 



have faded out. Most of the nonthermal emission is so smoothed by cosmic ray 
transport that the spatial distribution of its sources cannot be deduced in detail. 
Thus a various number of other possible sources for cosmic ray acc eleration has b een 
proposed, but none of them was really convincing (a review is in ICondonlFmoi . A 
five-band (radio, FIR, near infrared, blu e, and X-ray) study of luminosity correla- 
tions in spirals led Fabbiano et al. ( 1988| ) to conclude that the radio emission from 
starbursts must originate in the young stellar population but no clear conclusion 
can be drawn from the luminosity correlations in less active spiral galaxies. 



1.3 The radio/FIR correlation 



A correlation between the A = 10 u mid infrared and 14 15 MHz radio luminosities of 



Seyfert galaxies was discovered bv Ivan der Kruitl f|l97lh and soon extended to normal 



spiral galaxies (jvan der Kruiti 11973 ). At first both the infrared and radio emission 
were thought to be synchrotron radiation from relativistic electrons accelerated by 
nuclear monsters (e.g. massive black holes in Seyfert galaxies or other AGN). Then 



Harwit fc Pacinil (ll975h proposed that the infrared is thermal reradiation from dusty 



Hn regions, while the 1415 MHz is dominated by synchrotron radiation from rel- 
ativistic electrons accelerated in SNRs from the same population of massive stars 
that heat and ionize the Hn regions. This is still the current interpretation (cfr. 
Sect. Q and Q)) . 

The real significance of the FIR/radio correlation for normal galaxies — its tight- 
ness and universality — was not appreciated until the large IRA S survey appeared 
( Dickev fc Salpeterl Il984l be Jong et aD Il98,4 iHelou et all Il985l ) . The correlation 



has a linear slope and spans ~ 5 orders of magnitude in luminosity with less than 50% 



dispersion (|Yun et al.ll2001t iWunderlich et al.l 119871 ). making it one of the tightest 



objects, and holds up to z 




The correlation not only holds on 



global scales but is also found to h old within the disks of galaxies, down to scales of 
the order of a few hundreds of pc ( Beck fc Golla 1988| ). 

What is most extraordinary about this rela tionship is that its slope is unity (the 
most up-to-date result gives L\ A oc LgQ 9 ^ ±0,01 . mm et al.ll2001l ). although it couples 



four different radiative processes: i) thermal emission in the FIR from warm dust 
around Hn regions; ii) thermal emission in the FIR from cool dust (the 'cirrus') 
associated with dust heated by the old stellar population; iii) synchrotron emission 
at radio wavelengths; iv) thermal free-free emission at radio wavelengths. 

Several studies have thus been d evoted to investig ate the correlation by focusing 
on one or more of these processes. Cox et al. (|l988l ) correlated the luminosities at 



60/i and 151 MHz (where the contribution of the free-free emission to the radio 
luminosity should be negligible) and found a non-linear slope of L FIR oc L?m 6±0,0 3 



After deconvolving the FIR luminosity into two temperature components (|Fitt et al 
1988j), the best fits yielded -^FiR,warm oc L^ ±om f or the warm component, and 



£fir,coo1 oc . Thus they suggested that the warm dust component was correlated 
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1.3 The radio/FIR correlation 
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Figure 1.7: The radio/ FIR correlation fo r a sample of 1809 galaxies from the 
IRAS 2 Jy catalogue, from lYun et al. ( 2001 ). The correlation spans five orders of 
magnitude in luminosity with a best fit slope of 0.99 ± 0.01. 



with the nonthermal radio luminosity, while the link between the cool dust and the 
synchrotron radiation remained unexplained. 

However, decoupling the warm and cool components from IRAS data alone should 
be taken with care, because the cool dust becomes visible only at wavelengths longer 
than 100/x. The PHOT instrument onboard the ISO satellite was better suited for 



Pierini et all (|2003h thus 



this task, because of the availability of a 170/i channel, 
used PHOT data taken at 60, 100 and 170/x for a sample of 72 galaxies and found 
(note that here 'FIR' indicates the bolometric luminosity for the dust) 



L1.4 oc L 



1.10±0.03 
FIR 



L\ A oc L 



1.03±0.03 
FIR, warm 



Ll.4 OC L 



1.13±0.04 
FIR,cool 



(1.1) 



thus confirming the view that the warm dust and the nonthermal radio luminosity 
are linearly correlated. Also their b est- fit slope for th e cool dust /radio correlation 



is consistent, within the errors, with iFitt et (jl988h 



Theoretical models for the radio/FIR correlation 

One of the earlie s t theories, the 'optica lly-thick' or 'calorimeter' theory (Volk 1989; 
Volk fc XiJll993 iLisenfeld et ap Ejjj), is based on three main assumptions. First, 
that all the far ultraviolet (i.e. ~ 5-13.6 eV) radiation from massive stars is absorbed 
by the dust grains within a galaxy. Second, that the energetic electrons produced 
by the supernova explosions of these stars lose most of their energy within the 
galaxy due to synchrotron and inverse Compton process. As both these processes 
are proportional to the number of massive stars, these calorimetric assumptions lead 
to the linear correlation. Finally, the tightness of the correlation is provided by the 
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1 Light emission from star forming galaxies 



third assumption that the energy density of the interstellar radiation field, f/ ra d, is 
in a constant ratio with the magnetic field energy density, Ub- 



An alternative theory put forward by Helou fc Bicav ( 1993j ) assumes the opposite 
extreme, an 'optically-thin' model, in which the cosmic rays and UV photons both 
have high escape probabilities. To provide the correlation in these optically- and 
cosmic ray thin galaxies they rely on the assumptions that the UV or dust heating 
photons and the radio emitting cosmic rays are created in constant proportion to 
each other, which is again related to star formation. Then, to obtain a linear cor- 
relation, it is required by the theory that the magnetic field strength and the gas 
density are coupled. 



A challenge to both these theories is the model bv lNiklas fc Beckl ([1997D . In this 
work they argue that observations indicate that within most galaxies, the cosmic ray 
electrons lose very little energy before they escape. These same galaxies are optically 
thick to UV photons, thus both the c alorimetric and optica lly thin models are not 
supported by the observations. In the iNiklas fc Beck! (|l997l ) model, the controlling 
factor they put forward for the correlation is the gas density: they assume that both 
the star formation rate (and thus dust heating) and the magnetic field strength 
(which determines the synchrotron emission) are correlated with the gas volume 
density. 

All these models require a sort of fine tuning between the magnetic field strength 
and another parameter (the radiation field or the gas density) . As a way out from the 
fine tuning hypoteses, two mechanisms have been proposed. iBressan et al. ( 20021 ) 
suggested that, in a 'calorimeter' theory, the radio/FIR correlation may arise if the 
SFR remains contant for a time longer than the synchrotron-loss timescale. Under 
this hypothesis, the nonthermal radio luminosity of a galaxy is proportional to the 
integral of the synchrotron power over the electron lifetime, and an increase of the 
former in a larger magnetic field is compensated by a shortening of the latter. 



Ano ther me chanism was p r opose d bv lGroves et al. (2003). Following the lNiklas fc Beck! 
(1997) model, Groves et al. ( 2003| ) proposed magneto-hydrodynamic turbulence as 
a possible mechanism to provide coupling between the gas density and the magnetic 
field strength, under the requirement that the radio and FIR emission is produced 
in the same volume element. 



A 'definitive' model has not yet emerged. Further work is needed, as the current 
models do not explain the (non-linear) correlation between the radio emission and 
the cool dust emission. Also, the possibility that both the FIR and radio emis- 
sion increase more than linearly with increasing SFR has been invoked 
Possible reasons might be that low luminosity galaxies have significantly less dust 
absorption (hence their infrared emission misses most of the star formation), and 
that the thermal fraction of radio emission increases with the SFR. Thus the ra- 
dio/FIR correlation would arise as a sort of conspiracy. 



Be3|2003). 
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1.3 The radio/FIR correlation 
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Figure 1.8: A co mparison of S FR es timates from radio, infrared, Ha and ultra- 
violet emission, from Cram et al. ( 1998| ). The radio/FIR correlation (upper panel) 
has the least scatter. 
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1 Light emission from star forming galaxies 



1.4 SFR indicators 

Here we will present some very simple considerations, on which the models reviewed 
in Sect. 11.31 are based, on how to link the luminosity and the star formation rate. 
If we assume that the dust surrounding the SF region is optically thick to the UV 
radiation emitted by young stars, and there is no other process leading to far infrared 
emission, we may write 

L FIR oc/ kty(M)r(M)L(M)dM (1.2) 

J Mi 

where \I/(M) is the (time independent) Initial Mass Function (IMF), i.e. the proba- 
bility that a star of mass M is formed in an episode of SF, k is a normalization for the 
IMF giving the actual number of stars formed per unit time, and r(M) and L(M) 
are the lifetime and UV luminosity of a star with mass M, respectively. Under the 
hypothesis of continous star formation, or if a burst is long enough to reach a steady 
state in the mass distribution of stars, then k^f(M)r(M) is the number of stars with 
mass M actually present. M\ may be taken as the lowes t possible mass for a star 
which emits most of its luminosity in the UV band (6M Q , iDevereux fc Youngl ll990: 
however in the following the SFR calibrations will be given for M\ = 5M ). M u is 
upper bound of the IMF; it is usually assumed M u = 1OOM or 12OM , however 
for a Salpeter IMF (\I/(M) oc M~ 2,35 ) the integral is not really dependent on this 
parameter. The life time for a main sequence star may be written as 

qMX 

roc— — (1.3) 

where q is the hydrogen fraction actually burned, and X is the hydrogen abundance. 
Thus L(M)t(M) oc M. Since we define 



SFR(M>Mi)=/ Mkm(M)dM (1.4) 

J Mi 



then we get that L FIR ocSFR. 

Thus we may take ( Cram et al. 1998) 



SFR 60M (M > 5M ) = 2.0 ■ KT 31 L 60M (1.5) 
with Lqq^ in erg s _1 Hz -1 . Another calibration ( Kennicutt 19981 : Cram et al.l [l998) 



SFR FIR {M > 5M ) = 4.5 • 10~ 44 L Fm M yr" 1 (1.6) 
makes use of the FIR parameter! Helou et al. 19851 ) 



FIR = 1.26 • 10- 11 (2.58S' 60M + S 1Q0 ^ ergs" 1 cm" 2 (1.7) 
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with Seon and Sioo^ in Jy. The FIR parameter estimates the flux which would be 
measured by a theoretical detector with flat response between 42.5// and 122. 5/i if 
the underlying spectrum is a modified black body law 



B(u,T) 



2hu 3 



1 



hv 

ekr 



1 



1-e 



;i.s) 



with 20K < T < 80K and < f3 < 2. 

If, on the other hand, we assumed that the dust and gas are optically thin for the 
emission in the ultraviolet band and the Ha line, a similar reasoning might have been 
made by changing Lfir with Luv or L-R a in Eq. (jl.2|) . Balmer line emission from 



star- fo rming galaxies is formed in H II regions ionized by early-type stars. iKennicutt 
( 1983h has determined the theoretical relationship between the Ha luminosity and 
the current rate of star formation in a galaxy in a form corresponding to 



SFR Ha (M > 5M e 



QJ 



1.5 • 10 41 erg s" 



M yr" 



;i.9) 



The use of far-UV or U-band luminosities to infer the current SFR rests on the 
idea that UV emission contains a substantial contr ibution of light from the photo- 
spheres of young, massive stars ((Cowie et al.lll997l ). However, the calibration of a 
UV indicator is less straightforward than those listed above, because of the large 
numbers of relatively old stars that also c ontribute to t he act ual U-band luminosity 
of disk galaxies. To avoid this problem, ICram et al. I (ll99Sl) used the relationship 



between SFR and far-UV luminosity given by iCowie et al 



19971 ). scaled from 250 



nm to the U-band using relevant synthesized spectra, to derive 



SFR UV (M > 5M G 



©J 



1.5 • 10 29 erg s" 1 



M Q yr- 



;i.io) 



The two hypotheses made above (the dust and gas being optically thick or thin to 
UV radiation) are obviously at odds. Here we will not venture into details, but rather 
remind that in a patchy distribution of dust /gas clouds (a "leaky cloud" model, 
Devereux fc Young|l99Cll : Boulanger k, Peraultlll988| ). a fraction of the ultraviolet /Ha 
radiation produced in an Hn region may still be able to escape. Clearly, extinction 
corrections will be needed to recover the SFR from UV/Ha observations. SFR es- 
timates from UV and Ha vs. radio are shown in Fig. (jl.8|) . middle and lower panel 
respectively. If we assume that the radio luminosity is a good SFR indicator (see 
next paragraph), then the high scatter and the non-linearity of the (UV, Ha)-radio 
relation cleary shows the effects of extinction. Note that Eq. (11.91) incorporates a cor- 
rection for the average extinction suffered in spiral galaxies ( Kennicuttlll983j ). thus 
the scatter is entirely due to the individual variations of the amount of absorption 
in individual galaxies. 



Th e radio emission may be linked to the SFR with similar arguments ((Con don 
1992|. If all stars with mass M > M sn end their life as core-collapse supernovae, 
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then the supernova rate may be defined as 

«Af u 



fctf (M) dM 



A/ S N 



and, assuming that SNR are radio sources only during their adiabatic lifetime, the 
nonthermal luminosity might be written as 



T P- 1 /17„-2/17 1 . 



;i.i2) 



where i£ is the supernova explosion energy, and n is the ambient particle density 
However, it is impossible to reconcile the supernova rates needed by Eq. (J1.12j) with 
those obse rved in galax i es, since the former are abo ut ten times larger. Thus it was 
proposed ( Poolev 19691 : Ilovaiskv fc Lequeux 19721 ) that SNRs may accelerate the 
relativistic electrons producing the nonthermal radio emission, but > 90% of this 
emission must be produced long after the individual supernova remnants have faded 
out and the electrons have diffused throughout the galaxy Thus it may be better 
to use the empirical relation observed in the Milky Way 



;i.i3) 



whit Lnt in units of 10 2 ^ erg/s, y in G Hz and z^sn in yr _1 ; the average a found in 
spiral galaxies is ~ 0.75 ()Condonlll992i ). Eq. (jl.13)) probably applies to most normal 
galaxies, since significant variations in the ratio Z/nt/^sn would violate the observed 
radio/FIR correction. Also I vfilk et ni l fTiwiil argued that the cosmic ray energy 
production per supernova is the same in M82 as in the Galaxy. Thus we have, in 
terms of the SFR: 



5.3 ■ 10 2 V 



Y SFR. 



;i.l4) 



1.5 X-ray emission: the Chandra and XMM- Newton 
revolution 

It has been known since the late eighties that star forming galaxies are luminous 
sources of X-ray emission, due to a number of X-ray binaries, young supernova 
remnants, and hot plasma associated to star forming regions and galactic winds 
(jFabbianol Il989h . The X-ray emission is a key tool in understanding the star for- 
mation processes, since it traces the most energetic phenomena (mass accretion, 
supernova explosions, gas heating, cosmic rays acceleration), and it is a probe into 
the dense and dusty star forming regions. 

In the last few years — i.e. after the launch of Chandra and XMM-Newton — , 
the 'common picture' of normal galaxies underwent a dramatic review, thanks to 
the quantum leap in sensitivity, angular and spectral resolution achieved by the 
new observatories (cfr. Table ITTj) . As an example, in lPlate Ilia the images of the 
Antennae galaxy (NGC3038/9) taken with ROSAT and ASCA are shown beside an 
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Telescope 
(instrument) 


PSF 
(FWHM) 


Effective area Spectral resolution 
(cm" 2 at 1 keV) (E/AE at 1 keV) 


Bandpass 
(keV) 


Einstein (IPC) 


60" 


100 


0.7 


0.4-3.5 


ROSAT (PSPC-B) 


25" 


200 


2.5 


0.1-2.4 


ROSAT (HRI) 


3" 


80 




0.1-2.4 


ASCA* 


50" 


420 


10 


0.5-10 


BeppoSAXt 


90" 


60 


5 


0.1-300 


Chandra (ACIS) 


0.5" 


350 


10 


0.5-8 


XMM-Newton (EPIC)* 


6" 


1500 


13 


0.5-10 


XMM-Newton (RGS) 


(~ 15") 


100 


300 


0.3-2.0 



Table 1.1: Characteristics of past and current X-ray telescopes, * effective area: 
GIS+SIS; energy resolution: SIS. ^effective area: LECS; PSF: MECS. ^effective area: 
MOS+PN; energy resolution: MOS. 



optical image from the Digitized Sky Survey. The improvement provided by Chandra 
(jPlate IHb ) is striking, since it provides a ten-fold increase in angular resolution with 
respect to ROSAT, while still offering a spectral resolution comparable to ASCA. 

Point sources (X-ray binaries and young SNRs) have emerged as the dominant 
contributors to the overall X-ray emission, especially in the hard band (i.e. 2.0-10 
keV); diffuse emission is nonetheless significant in the soft band (0.5-2 .0 keV) for 
the galaxies with the largest SFRs, such as M82 ( Gri ffiths et al. 2000) where the 



diffuse emission accounts for ~ 30% of the soft emission. In the following, the main 
characteristics of the different kinds of sources will be briefly discussed. 

X-ray bi naries are the brightest class of point sources in galaxies (rev iews may 
be found in White et al. 19951: Ivan Paradiisl 19981 : IPersic k, Rephaeli 20021 ) ; they are 
formed by a collapsed star (black hole or neutron star) with a normal companion 
(main sequence or giant). Although the properties of their X-ray emission (lumi- 
nosity, variability, spectrum, whether the emission is persistent or transient) depend 
on a number of parameters (mainly the masses of the two stars and the accretion 
rate), an essentially full characterization of X-ray binaries can be made on the basis 
of the mass of the donor star. We have High Mass X-ray Binaries (HMXB) when 
the donor is a main sequence star with M > 8M Q , and Low Mass X-ray Binaries 
(LMXB) when the donor is a post main sequence star with M < 1M Q (binary 
systems with intermediate mass donor stars turn out not to be effective X-ray emit- 
ters). This main distinction has to do with the nature of the mass transer — wind 
accretion in HMXB, transfer via Roche-lobe overflow in LMXB. 

The spectral shape of X-ray binaries may be basically seen as the sum of two 
components: a black body with kT > 1 keV ('soft component', associated with an 
accretion disc) and a power law S(E) oc E l ~ T with slope T ~ 1.2 for the HMXB and 
~ 1.4 for the LMXB ('hard component', arising from Comptonization of the photons 
form the accretion disc); the relative weight of the two components may vary with 
time, and one of the components may even completely dominate the overall emission 
- a transition occurs between the so-called soft (or high) and hard (or low) states 
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(cfr. Sect. I2.fij) . The overall X-ray luminosity is usually comprised in the 10 37 -10 39 
erg/s interval, with the HMXBs somewhat brighter on average, and depends for 
both types on how close the normal star is to fill its Roche lobe: the closer, the 
brighter. 

Supernova remnants are moderately luminous sources (Lx ~ 10 37 erg/s) in the 
early part of their life, namely during the free expansion phas e and the beginning of 
the adiabatic expansion phase (|Woltierlll972l IChevalierllllmh . Their spectrum can 
be described as a thermal plasma with decreasing temperature, from kT ~ 40 to 
~ 0.4 keV in about 10 3 yr, as the remnant expands (assuming a rather dense ISM 
as should be th e case in the most actively star for ming galaxies; namely n ~ 10 3 
cm -3 , following Chevalier Sz Fransson 2001 : see also lDraine &: Woodsl 1991 ). This 
rather short lifetime makes SNRs a minor component among X-ra y sources in galax- 
ies, accounting for < 10 % of the total number of point sources ([Pence et al.l 12001: 
Persic fc Rephaelill2002h . It is nonetheless interesting to note that since the radio- 
emitting phase only begins 10 4 yr after the supernova explosion, X-ray and radio 
SNRs represent two different evolutive phases of the same population of objects; in- 
deed the comparison of Cha ndra and radio images has sh own that sources emitting 
in both bands are very rare (|Chevalier Sz FranssonlEoOlh . 

The interstellar medium (ISM) itself, when heated by supernova explosions, turns 
into an X-ray source. Soft, diffuse emission in galaxies was already observed with 
Einstein, and at least one thermal component with a temperature around ~ 1 ke V 



was r equired to fit the ASCA and BeppoSAX data fPt ak et al.lll997t ICappi et al. 



1999). Its spectrum may be described as a Bremsstrahlung component plus a number 



of atomic recombination lines. Depending on the ISM density, the SFR, and the 
gravity potential of the galaxy, the hot high-pressure plasma may eventually find a 
way out of the galaxy disc and escape into the intergalactic medium — thus creating 
an outflow (jPlate V|) . Observations conducted with the ROSAT HRI instrument 
showed a spatial overlap of t he X-ray emission with Ha emission from the outflow 
f De lia Ceca et al.lll996lll997l ). It was not clear, however, if the X-ray emission came 
from the interface between the ambient material and the wind, or from the wind 
itself. Chandra ob servations showed tha t the X-ray emission has a conical, shell-like 
shape (IPlate Vlb : fstrickland etail2 000). It can be explained as the hot rarefied and 
expanding wind heats the surrounding material which in turn produces the X-ray 
emission. XMM- Newton observations (jPlate Vlh ) have shown that also the galactic 
disc is an X-ray emitter. An interesting and yet unresolved question is the fate of 
the metals produced in the supernova explosions: whether they mix with the ISM 
or escape from the galaxy; cfr. Chap. O 

The total spectrum of a star forming galaxy is thus made up by the sum of a 
thermal spectrum from hot plasma and of a power law representing the emission from 
the X-ray binaries. While a detailed discussion of the relative contribution of the 
differe nt classes of sources to the total spectrum may be found in lPersic fc Rephaeli 
we will focus here on two simple considerations: i) the power law component 
is often found to be absorbed, with column densities aroud io 21 ~ 22 cm -2 ( Ptak et al 



1999; iMoran et Il999l ^ : ii) the thermal component has a steeper slope than the 
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channel energy (keV) 
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Energy 



Figure 1.9: X-ray spectrum of M82. Upper panel: observed data from the pn 
instrument onboard XMM-Newton, rebinned to achieve a significance of at least 15a 
per energy bin. The continous line shows the model, folded with the instrumental 
response. Lower panel: model spectrum (continous line) composed by a hot plasma 
(dashed line) plus an absorbed power law (dot-dashed line). 
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power law. The joint effect of the steeper slope of the thermal component, and of 
absorption on the power law, makes the two components mainly emit in different 
energy bands: the plasma dominates at softer energies (< 1-2 keV) and the binaries 
at harder ones. As an example, the spectrum of the galaxy M82 obtained with the 
pn instrument onboard XMM- Newton is shown in Fig. It is interesting to 

notice that, although the individual spectrum of X-ray binaries has a slope around 
r ~ 1.2-1.4, the resulting total spectrum of a galaxy may have a flatter slope if 
a number of binaries is present with different absorbing column densities; cfr. the 
Chandra spectrum of the point sources in M82, Fig. 15.31 Conversely, if the spectrum 
of the X-ray binaries is not heavily absorbed, the more intense the plasma emission, 
the more steep the overall spectrum of a galaxy. As a consequence, a relatively 
large spread may be found in the average hard X-ray spectrum of starburst galaxies 
( Ptak et al.lll999h . 



An active galactic nucleus (AGN) may also be present in star forming galaxies. 
While it has been proposed that nuclear activity might be related to enhanced star 
formation (the 'starburst /AGN connection') as a consequence of the disposal of large 
quantities of gas, nuclear activity is not truly a star formation related process. If 
present, moreover, an AGN tends to dominate over the X-ray emission from the rest 
of the galaxy. A typical AGN luminosity (L x > 10 42 erg/s) exceeds by (at least) 
a few orders of magnitude the emission from a normal galaxy (Lx ~ 10 40±1 erg/s). 
Thus the main concern about AGN in X-ray studies of starburst galaxies has been 
so far to exclude them from the galaxy samples (cfr. Sect. 12. 1|) . The high resolution 
of Chandra allows to resolve most nearby g alaxies, so that AGN may be identified 
and properly accounted for in data analysis ()Ho et al.ll200l[ ). However, the problem 
still arises for more distant, hence unresolved, galaxies. 

Optical spectroscopy is the main tool used in the classification of Seyfert vs. star 
forming galaxies. The equivalent widths of emission lines from the centre of a galaxy 
are used in diagnostic diagrams to determine the excitation and ionization states of 
the nebular gas, which in tu rn are determined by the nature of the energy source 
(jVeilleux OsterhroclJ Pi 987t ). Although a determination of the X-ray luminosity 



and/or an optical spectrum usually suffice, this can prove to be a difficult task, 
especially for luminous dusty objects (such as NGC 3256, NGC 3690, NGC 6240, 
Arp220) in which the overall luminosity is still sustainable by star formation, and 
AGN signatures may be hidden by a bsorption. NGC 3690 ( Delia Ceca et al. 2002h 



and NGC 6240 ( Vignati et al. 19991 ) are emblematic cases: both of them are cur 



rently undergoing a powerful starburst (100-200 M /yr), and the optical spectrum 
does not show clear AGN signatures; however when observed with the PDS instru- 
ment onboard BeppoSAX they have shown the emergence at energies larger than 
10 keV of a heavily absorbed (iVn ~ 10 24 cm -2 ) nuclear component. 

At high redshifts (z ~ 1), the use the diagnostic diagrams in separating Seyfert 
galaxies from star forming ones has to be treated with caution because of two ef- 
fects: a possibly lower signal/noise ratio, and the shift of emission lines outside the 
observing spectral band. Moreover, because of the little angular size (a few arcsec) 
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of high redshift galaxies, it is not possible to obtain a spectrum of the central regions 
only. Thus the signatures of a low luminosity AGN may be diluted in the overall 
radiation of the galaxy. With this caveat, optical spectroscopy may still be regarded 
as an effective method in selecting candidate star forming galaxies at high redshifts. 
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2 The X-ray luminosity as a SFR 
indicator 



While radio continuum and far infrared (FIR) luminosities of star-forming galaxies 
are known to show a tight linear relationship spanning four orders of magnitude in 
luminosity and up to a redshift ~ 1.3 (Sect. 11.3)1 . a relation between FIR and X- 
ray luminosities was found in the early nineties but its details remained somewhat 
controversial. A non linear (Lx oc ^fir) an< ^ highly scattered (dispersion of about 
2 dex) relation was found between FIR and soft (0.5-3.0 keV) X-ray luminosities of 
IRAS-bright and/or interac ting/peculiar galaxies measured by the Eins tein satellite 
( Griffiths fc Padovani 1990l ) . A somewhat different result was found by David et al. 



(1992J), i.e. a linear relation between FIR and 0.5-4.5 keV luminosities for a sample 
of starburst galaxies observed by Einstein. A la rge number of upper limits to the 
X-ray flux (12 upper limits vs. 11 detections for Griffi ths fc Padovani 1990f ) along 
with high uncertainties in the X-ray and FIR fluxes may explain this discrepancy. 
Moreover, these studies suffered by the lack of knowledge about spectral shapes and 
internal absorption in star forming galaxies caused by the limited sensitivity and 
spectral capabilities of the IPC detector onboard Einstein. 

In this chapter, with the high sensitivity and the broad-band spectral capabilities 
of the ASCA and BeppoSAX satellites, we extend these studies to the 2-10 keV band 
which is essentially free from absorption. In the following paragraphs a sample of 
nearby star forming galaxies is assembled (Sect. 12. 1J) and linear relations among 
radio, FIR and both soft and hard X-ray luminosities are found fSect. l2~2~j) . Possible 
biases are discussed and the use of X-ray luminosities as a SFR indicator is proposed 
(Sect. 12.3)1 . In Chap. El we present a study of star-forming galaxies in the Hubble 
Deep Field North and test the validity of the X-ray SFR law. Implications for 
the contribution of star-forming galaxies to the X-ray counts and background are 
discussed in Chap. |U 

Throughout this thesis we assume, unless otherwise stated, H = 70, Qm = 0.3 
and = 0.7. 



2.1 The local sample 

The atlas of optical nuclear spectra by Ho et al. ( 19971 ) (hereafter HFS97) represents 
a complete spectroscopic survey of galaxies in the R evised Shapley-Ames Catalog 
of Bright Galaxies fRSA; ISandage fe Tammann 1981) and in the Second Reference 
Catalogue of bright galaxies (RC2; de Vaucouleurs et al. 19761 ) with declination 5 > 
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0° and magnitude Bt < 12.5. O ptical spectra are classified in HFS97 on the basis of 



j ptic 

line intensity ratios according to lVeilleux fc Osterbrockl (|l987f ): galaxies with nuclear 



line ratios typical of star-forming systems are labeled as "H II nuclei" . This sample 
of Hn galaxies contains only spirals and irregulars from Sa to later types, except 
for a few SO which were excluded from our analysis since their properties resemble 
more those of elliptical galaxies. 

A cross-correlation of the HFS97 sample with the ASCA archive gives 18 galaxies 
clearly detected in the 2-10 keV band with the GIS instruments. Since most of these 
galaxies were observed in pointed observations (rather than in flux-limited surveys) 
it is unlikely that the results hereafter shown are affected by a Malmquist bias. Four 
additional objects in the field of view of ASCA observations were not detected: the 
2-10 keV flux upper limits are too loose to add any significant information, and thus 
we did not include them in the sample. The cross-correlation of the HFS97 sample 
with the BeppoSAX archive does not increase the number of detections. When a 
galaxy was observed by both satellites, the observation with better quality data was 
chosen. 

Far infrared fluxes a t 60m and lOOu wer e taken from the IRAS Revised Bright 
Galaxy Sample fRBGS. ISanders et al.l2003f ) which is a reprocessing of the final IRAS 
archive with the latest calibrations. While the RBGS measurements should be more 
accurate, we checked th at the use of the older catalogue of IRAS observations of 
large o , coupled with the Faint Source Catalogue 

(FSC, iMoshir et al.l Il989l ) for smaller galaxies, does not signi ficantly chan ge our 



statistical analysis. FIR fluxes for NGC444 9 were taken from Rush e t al. 



Radio (1.4 GHz) fluxes were obtain ed from the Condon et al. ( 1990Lll996| ) catalog 



mes 



(excep t for NGC4449, taken from Havnes et al.ll975h . Distances were taken from 
and corrected for the adopted cosmology. 



Part of the X-ray data have already been published; in the cases where published 
data were not available in a form suitable for analysis, the original data were re- 
trieved from the archive and reduced following standard procedures and with the 
latest available calibrations. Images and spectra were extracted from the pipeline- 
screened event files. The images were checked against optical (Digital Sky Survey) 
and, where available, radio (1.4 GHz) images in order to look for possible source 
confusions. Fluxes were calculated in the 0.5-2.0 and 2-10 keV bands from best-fit 
spectra for the GIS2 and GIS3 instruments and corrected for Galactic absorption 
only. The uncertainty on the fluxes is of the order of 10%. Depending on the quality 
of data, the best-fit spectrum is usually represented by a two-component model with 
a thermal plasma plus a power-law or just a power-law. 

The galaxies IC 342 and M82 have shown some variability, mainly due to ultralu- 
minous X-ray binaries. For the two of them we summarize in Sect. 12.61 the results 
from several X-ray observations and estimate time-averaged luminosities. 

One object (M33) was not included in the sample since its broad-band (0.5-10 
keV) X-ray nuclear spectrum is dominated by a strong variable source (M33 X-8) 



1 Blue-light isophotal major diameter (-D25) greater than 
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identified as a black hole candidate (jParmar et alJEoOl ). 

Therefore, the sample (hereafter local sample) consists of the 17 galaxies listed in 
Table 12.11 Since it is not complete in a strict sense due to the X-ray selection, we 
have checked for its representativeness with reference to the SFR. The median SFR 
value for HFS97 was computed from the FIR luminosities and resulted SFR mo d ~ 
1.65 M /yr. Considering objects with SFR > SFR med there are 14 galaxies in the 
local sample out of 98 in HFS97 (14%), while there are 3 objects with SFR < SFR me d 
(3%). Thus the high luminosity tail is better sampled than the low luminosity one. 

We also include data for 6 other well-known starburst galaxies which were not in 
the HFS97 survey because they are in the southern emisphere. On the basis of their 
line intensity ratios 2 they should be classified as H II nuclei. In Table 12.11 we label 
them as supplementary sample. 



2.2 The radio/FIR/X-rays correlation 



As a preliminary test, we perform a least-squares analysis for the well-known ra- 
dio/FIR correlation, which yields 



Log(L FIR ) = (0.98 ± 0.06) Log(Li. 4 ) + 15.4 ± 1.6 



(2.1) 



The dispersion around the best-fit relation is given as the estimate s of the stan- 
dard deviation a: 



1 



N — v 



^2 ( L °g ^obs - Log L 

pred ) 



(2.2) 



where v is the number of free parameters and N is the number of points in the fit), 
L pred is the luminosity expected from the best fit relation and L obs the observed one. 
For the radi o /FIR correlation ( Eg. I2.1J) we find s = 0.18. 

Following iHelou et al.l (|l985T ) we also calculate the mean ratio q between the 
logarithms of FIR and radio fluxes, obtaining q ~ 2.2 ± 0.2. This value is consistent 
with the mean q = 2.34 ± 0.01 for the 1809 galaxies in the IRAS 2 Jy sample by 
Yun et all f)200lh . 



Soft X-rays 

A test for the soft X-ray/FIR/radio relations fFig. 12. lj) yields 



Log(L .5- 2 ) = (0.87 ± 0.08) Log(L FIR ) + 2.0 ±3.7 (2.3) 
Log(L . 5 „ 2 ) = (0.88 ± 0.08) Log(L L4 ) ± 14.6 ±2.2 (2.4) 

References: NGC55 ^Webster fc Smdthl Jl98.^: NGC253 1672 fc 1808 -iKewlev et all (l200lh : 
NGC 3256 - iMoran et all (jl999^ : Antennae - iRubin et alJ l|l970|) , iDobrodiv fc Pronikl 1)1979^ 
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Fluxes and Luminosities: Main Sample 



Galaxy 


D 


Fo.5-2 


£o.5-2 


-P2-10 


£2-10 


Ffir 


Lfir 


Fia 


■£'1.4 


Refs. 


M82* 


5.6 


97 


3.6 


290 


11 


67 


25 


7.7 


2.9 


1 


M101 


5.7 


5.4 


0.22 


6.8 


0.27 


6.0 


2.4 


0.75 


0.30 


this work 


M108 


15 


4.4 


1.2 


6.0 


1.6 


2.0 


5.4 


0.31 


0.83 


this work 


NGC891 


10 


8.3 


0.99 


19 


2.3 


4.5 


5.4 


0.70 


0.84 


this work 


NGC1569 


1.7 


5.4 


0.019 


2.2 


0.0077 


2.5 


0.088 


0.41 


0.014 


2 


NGC2146 


18 


8.2 


3.4 


11 


4.5 


7.3 


30 


1.1 


4.5 


3 


NGC2276 


39 


2.1 


3.9 


4.4 


8.1 


0.85 


16 


0.28 


5.2 


this work 


NGC2403 


4.5 


16 


0.39 


9.3 


0.23 


2.7 


0.65 


0.33 


0.080 


this work 


NGC2903 


6.7 


7.9 


0.43 


7.0 


0.38 


3.7 


2.0 


0.41 


0.22 


this work 


NGC3310 


20 


7.4 


3.5 


2.1 


1.0 


1.7 


8.1 


0.38 


1.8 


4 


NGC3367 


46 


1.8 


4.5 


1.6 


4.0 


0.38 


9.5 


0.10 


2.5 


this work 


NGC3690 


49 


5.7 


17 


11 


32 


5.3 


150 


0.66 


19 


4 


NGC4449 


3.2 


8.3 


0.10 


4.8 


0.060 


1.9 


0.23 


0.6 


0.074 


5 


NGC4631 


7.1 


9.4 


0.57 


9.3 


0.57 


4.9 


3.09 


1.2 


0.73 


6 


NGC4654 


18 


0.6 


0.2 


0.9 


0.3 


0.93 


3.5 


0.12 


0.46 


this work 


NGC6946 


5.9 


30 


1.2 


12 


0.49 


7.9 


3.2 


1.4 


0.57 


this work 


IC342 


4.2 


18 


0.38 


110 


2.3 


11 


2.3 


2.3 


0.49 


this work 



Supplementary Sample 



NGC55 


1.4 


18 


0.040 


6.8 


0.015 


4.7 


0.10 


0.38 


0.0084 


6 


NGC253* 


3.2 


25 


0.31 


50 


0.62 


49 


6.1 


5.6 


0.69 


1 


NGC1672 


16 


5.8 


1.7 


6.1 


1.8 


2.3 


6.8 


0.45 


1.3 


this work 


NGC1808 


11 


6.5 


1.0 


7.6 


1.2 


5.3 


8.2 


0.52 


0.81 


this work 


NGC3256 


40 


9.0 


17 


6.2 


12 


4.8 


92 


0.66 


13 


7 


Antennae 


27 


7.2 


6.3 


5.3 


4.7 


2.6 


23 


0.57 


5.0 


8 



Table 2.1: Data for galaxies in our local samples. All galaxies were observed 
with ASCA, except those marked with * observed by BeppoSAX. Distances in Mpc; 
X-ray fluxes in 1CT 13 
in Jy; X-ray luminosities in 10 
luminosities in 10 29 erg s _1 Hz 
the order of 10%. 



erg s cm , FIR fluxes in 10 erg s cm and radio fluxes 
erg s _1 , FIR luminosities in 10 43 erg s _1 and radio 
1 . The uncertainty on fluxes and luminosities is of 



References: 1 CaDDi et al. (1999T): 2 Delia Ceca et al. 


(Il999h: 4 Zezas et al. 


1998h: 5 Delia Ceca et al.1 (|l997) 


Moran et all (199flh: 8 


Sansom et al.l ( 1996f). 



( 19961): 3 



6 Dahlem et al. 



Delia Ceca et al. 
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2.2 The radio /FIR/X-rays correlation 



with s ~ 0.26 and 0.24 respectively. 

Our result is consistent with the L 5 _ 4 5 oc L FI 9 ^ ±0 ' 06 relation found bv lDavid et al. 



(1992) for normal and starburst galaxies from the IRAS Bright Galaxy Sample, 
but it is only marginal ly consistent with the much flatter and more dispersed re- 
lationship obtained by Griffiths fc Padovani ( 1990h for a sample of IRAS selected 
galaxies (£o.5-3keV oc Lg 6 ^ ±0 ' 14 ) and for a sample of starburst /interacting galaxies 

it „ r0.70±0.12>> 

Oo.5-3keV OC L mfl ). 

The inclusion of the objects of the supplementary sample (Table 12.1(1 does not 
significantly change the slopes, i.e. L 5 _ 2 . okev oc Lp^ ±0 ' 07 ; likewise, if we use the 
60/z luminosity instead of FIR, we obtain L0.5-2.0keV oc Lqq^ ±0 m . 

By assuming an exactly linear slope, the best fit relations for the local (lo- 
cal+supplementary) sample become: 



Log(L . 5 - 2 ) = Log(L FIR ) - 3.68 (3.70) (2.5) 
Log(L . 5 - 2 ) = Log(L!. 4 ) + 11.08 (11.10) (2.6) 

with s ~ 0.27 and 0.24 respectively. 

By applying the F-test we find that the free-slope fits are not significantly better 
than those with the linear slope, the improvement being significant only at the la 
level. 

Hard X-rays 

In Fig. (|2.2(1 we plot 2-10 keV luminosities versus FIR and radio ones. Least-squares 
fits yield: 

Log(L 2 _ 10 ) = (1.08 ± 0.09) Log(L Pm ) - 7.1 ± 4.2 (2.7) 
Log(L 2 _ 10 ) = (1.08 ± 0.09) Log(L L4 ) +8.8 ±2.7 (2.8) 

with s ~ 0.30 and 0.29 respectively. The linearity and the dispersion are not sig- 
nificantly changed neither by the inclusion of the supplementary sample (L2-10 oc 
L^ ±a07 and L 2 _ 10 oc L};° 1±007 , s ~ 0.27 and 0.26 respectively), nor by the use of 
the 6O/1 luminosity (L 2 _i oc L^ ±om ). 

By assuming an exactly linear slope, the best fit relations for the local (lo- 
cal+supplementary) sample become: 

Log(L 2 _ 10 ) = Log(L FIR ) - 3.62 (3.68) (2.9) 
Log(L 2 _ 10 ) = Log(Li. 4 ) + 11.13 (11.12) (2.10) 

with s ~ 0.29 for both fits. There is no significant improvement (less than la) in 
the free-slope fits with respect to the linear slope ones. 
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2 The X-ray luminosity as a SFR indicator 



SFR FIR M yr-> 

1 10 1 a 




10 41 10 4S 10 43 1 44 

L (KIR) erg s" 1 Hz~ 



10 45 10 46 




10 s7 10 aB 10 39 10 30 

L (1.4GHz) erg s 1 Hz 1 



10 3 



Figure 2.1: The 0.5-2.0 keV luminosity of local star-forming galaxies vs. radio 
and FIR ones. Squares: local sample; crosses: supplementary sample; dotted lines: 
eas. ff231l2~T)J) . 
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2.2 The radio/FIR/X-rays correlation 
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Figure 2.2: The 2.0-10 keV luminosity of local star-forming galaxies vs. radio 
and FIR ones. Symbols as in Fig. (J2HJ); dotted lines: Eqs. lj2.HI2.10j) . 
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2 The X-ray luminosity as a SFR indicator 



2.3 X-rays and the Star Formation Rate 

The existence of a tight linear relation implies that the three considered bands 
all carry the same information. Since the radio and far infrared luminosities are 
indicators of the SFR, the 0.5-2 keV and 2-10 keV luminosities should also be SFR 
indicators. However, before attempting to calibrate such relationships, we should 
c onsider the possi ble existence of selection effects. 

Ho et al. (1995) made a special effort in obtaining nuclear nebular spectra, so 



that a reliable spectral classification of the central engine could be derived. The 
main concern is the possibility that the Hn galaxies in the HFS97 sample could 
host a Low Luminosity AGN (LLAGN), which mig ht significantly contrib ute to the 
overall energy output. To check for this possibility, Ul vest ad &: Hoi ( 20021 ) observed 



with the VLA at 1.4 GHz a complete sample of 40 Sc galaxies in HFS97 with Hn 
spectra and did not find any compact luminous radio core. Instead, they found 
that the radio powers and morphologies are consistent with star formation processes 
rather than by accretion onto massive black holes; thus they suggest that H II nuclei 
intrinsically lack AGN. Therefore we believe that the HFS97 classification is reliable 
and that our sample is not polluted by AGN. 

It is also worth noticing that the soft X-rays relationships may involve some further 
uncertainties related to the possible presence of intrinsic absorption (negligible in the 
2-10 keV band for column densities usually found in normal galaxies). An example 
of this effect is the southern nucleus of NGC3256 (Sect. I2.5jh a dusty luminous 
merger remnant with two bright radio-IR cores where star formation is ongoing: 
while both of them fall on the radio/hard X-ray relation, only the northern core is 
on the radio/soft X-ray relation because the southern one lies behind a dust lane 
which absorbes at all wavelengths from ~ lu. to ~ 2 keV. The quasi-linearity of 
the soft X-ray relations suggests that absorption is unlikely to be relevant for the 
majority of the objects in our sample; however this effect may become significant at 
cosmological distances (z > 1 — 2) where galaxies have more dust and gas at their 
disposal to form stars. 

Thus we feel confident to propose the use of X-ray luminosities as SFR indica- 
tors. From eqs. (|2. 512.612. 912. HH) and from the calibrations of SFR^/^ and SFR radio 
(Sect. II. 4j) we derive: 

SFR = 2.2 • 10- 40 L . 5 _2kev M /yr (2.11) 
SFR = 2.0 ■ 10~ 40 L 2 _ 10kc v M /yr. (2.12) 

We also notice that there is growing evidence that star formation could play a 
major role even among those obj ects classified as LLAGN. The preliminary results 



of the Chandra LLAGN survey (|Ho et al.l 120 01) show that only about one third 



of LLAGN have a compact nucleus dominating the X-ray emission, while in the 
remaining objects off-nuclear sources and diffuse emission significantly contribute to 
the overall emission. 

Following this investigation, we have analyzed the relations betw e en ra dio /FIR/X- 
ray luminosities for the spiral galaxies in the iTerashima et al. (2002) sample of 
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2.3 X-rays and the Star Formation Rate 
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Figure 2.3: 2.0-10 keV vs. FIR luminosities of star-forming galaxies and LLAGN. 
Filled triangles: LINERs, filled circles: Seyfert galaxies. Open symbols as in 
Fig. (EID; line: Eq. flUJ). 



LLAGN, drawn from HFS97 and observed with ASCA, comprising 7 LINERs and 
15 Seyfert's with 4 • 10 39 < L 2 _ 10 < 5 • 10 41 erg s" 1 . We find that the X-ray/FIR 
and X-ray/radio luminosity ratios generally exceed those of star- forming galaxies, 
but about one third of the objects have ratios falling on the same locus of the 
star- forming galaxies (Fig. 12. 3|) . Therefore, the nuclear X-ray emission of these last 
LLAGN must be comparable to or weaker than the emission from star formation 
related processes. Moreover, the infrared (IRAS band) colours of these objects are 
also similar to those of star-forming galaxies, and completely different from those 
of QSOs, thus suggesting that the FIR luminosities of LLAGN may be powered by 
star formation. 
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2 The X-ray luminosity as a SFR indicator 



2.4 Comparison with other X-ray based SFR 
indicators 

Our inference of using the 2-10 keV lu minosity as a SFR in dicator is consistent with 
a st udy on Lyman-bre ak galaxies by Nandra et al. ( 20021 ) who have extrapolated 
the lPavid et al. (jl992l ) FIR/soft X-ray relation to the hard X-ray band obtaining a 
SFR/2-10 keV luminosity relation within 10% of our Eq. (|2.12|) . From a stacking 
analysis of Chandra data for a sample of optically selected Lyman-break and Balmer- 
break galaxies in the HDFN they find a good agreement of the average SFR as 
estimated from X-ray and extinction-corrected UV luminosities. 



A s imilar derivation of an Lx-SFR relation has been performed by Ba uer et al. 



(2002). using 0.5-8.0 keV and 1.4 GHz lum inosities of a joint sa mple of 102 nearby 



late-type galaxies observed with Einstein (jShaplev et al.l 120011 ) and of 20 galaxies 
w ith emission li n e spe ctra in the Chandra Deep Field N orth. 

Grimm et all (|2003l ): (see also iGilfanov et al.l l2004bh ave recently shown that 



the luminosity function of High Mass X-ray Binaries (HMXB) can be derived from 
a universal luminosity function whose normalization is proportional to the SFR. 
They also show that, due to small numbers statistics, a regime exists, in which 
the total X-ray luminosity grows non-linearly with the number of point sources. 
This has the consequence that also the total X-ray luminosity vs. SFR relationship 
is non-linear in the same regime, namely for SFR below 4 — 5M Q /yr. Although 
this might se e m at odds with our findings of Sect. I2.2| it is not really so, since 
Grimm et al. _f 2 QQ ^ |) analysis only refers to the contribution from HMXB. It has 
been shown ( Gilfan ov et al . 2004a) that once the Low Mass X-ray Binaries are taken 
into account, their result are in agreement with the X-ray/radio/FIR correlations 
we found. However, the number of LMXB should be related with the integrated star 
formation over a much longer time than the HMXB (e.g. several Gyr) and thus not 
trace the current SFR. Mantaining both the radio/FIR/X-ray correlations, and the 
interpretation of the emission in these bands as powered by recent star formation, 
should require that also the FIR and radio emission ar e non-line arly related with 
the SFR; this has indeed been suggested for the radio ([Belli l2003h . As to the FIR 
emission, further investigations will be required to clarify its relation with the warm 
and dust component; thus at present the question stays open. 



2.5 NGC 3256: a case for intrinsic absorption 



We present the test case of NGC 3256, a luminous dusty merger remnant in- 
cluded in the supplem enta ry sample. Detailed stu dies at s everal wavelengths (radio: 
Norris fc Forbeslll995l. IR: iKotilainen et al.lll996[ optical: iLmari et al1l2000L X-ray: 
Moran et al.l 19991 Lira et al. 19991 ) have shown that the energetic output of this 
galaxy is powered by star formati on occurring at s e veral locations, but mainly in 
the two radio cores discovered by Norris &: Forbes! ( 19951 ) and also detected with 
Chandra (IPlate Vllj) . 
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2.5 NGC3256: for intrinsic absorption 




Figure 2.4: A sequence of infrar e d ima ges of the centre of NGC3256 at increasing 



wavelenghts, from Kotilainen et al. (jl99fi| ). The southern nucleus lies behind a dust 



lane, and becomes visible in the L' band. The two nuclei are the brightest sources 
in NGC3256 in the radio and 2.0-10 keV bands (cfr. IPlate Vlljl . 
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2 The X-ray luminosity as a SFR indicator 



Author 




SFR ESTIMATE 


Method 


Ranalli et al. 


(2003J 


SFR= 
SFR= 


2.0 ■ 10- 40 L 2 _ 10 
2.2 ■ 10- 40 L . 5 _ 2 


radio/FIR/X-ray correlations, 
23 galaxies obs. w. ASCA 


Nandra et al. 


(2002) 


SFR= 


1.8 ■ 10- 40 L 2 -i 


0.5-4.5 keV/FIR correlation 
in David et al. (1992) 


Bauer et al. ( 




SFR= 


(1.7 ±0.5) • 10~ 43 - 

T 1.07±0.08 
'-^0.5-8 


102 galaxies with Einstein 
+20 galaxies in the CDFN 


Grimm et al. 


(2003J 


SFR= 
SFR= 


1.5 ■ 10- 40 L 2 _ 10 
(M > 4.5M /yr) 

(3.8 • 10- 40 L 2 _ 10 ) - 6 
(M < 4.5M /yr) 


theory + 23 nearby 
galaxies (mainly Chandra 
data) 



Table 2.2: A comparison of different SFR indicators based on the X-ray lumi- 
nosity. 



SFR L4 CH2 M s yr- 
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Figure 2.5: 2.0-10 keV vs. radio luminosities of the two nuclei of NGC 3256 
(filled circles). The upper and lower circles represent the north and south nucleus, 
respectively. The other symbols (same as in Fig. 12. lj) represent the other galaxies 
in the sample. 
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2.6 Effect of variability on the flux estimate 



The 3 and 6 cm radio maps ( Norris fc ForbesI 19951 ) reveal two distinct , resolved 
(FWHM ~ 1.2") nuclei and some fainter diffuse radio emission. Separated by 5" 
in declination, the two cores dominate the radio emission, the northern one being 
slightly (15%) brighter. They share the same spectral index (a ~ 0.8). Chandra 
observations (jPlate VII)) have shown that they also have similar 2-10 keV fluxes. 
Both of them follow the radio/hard X-ray correlation (Fig. 12. 5 j) . while only the 
northern one follows the radio/soft X-ray correlation. At other wavelengths the 
northern core is the brightest source in NGC 3256, while the southern one lies behind 
a dust lane and is only detected in the far infrared (A > lO um), as clearly shown i n 
the sequence of infrared images at increasing wavelengths in Kotilainen et al.l (|l996|) . 

Although the southern core appears as a bright source in the hard X-rays (E > 2 
keV), there are not enough counts to allow an accurate spectral fitting. However, it 
is still possible to constrain the absorbing column density by assuming a template 
spectrum, such as a simple power-law or the spectrum of the northern core, leading 
(after standard processing of the Chandra archival observation of NGC 3256) to an 
intrinsic N R = (2.2 ± 1.0) • 1 22 cm" 2 (JV H rc»i = 9.5-1 20 cm" 2 ), fully consistent 3 with 
the Ay = 10.7 estimated by Kotilainen et al. ( 19961 ) from infrared observations. We 
note that 84% of the 2-10 keV flux and only 10% of the 0.5-2 keV one are transmitted 
through this column density. Thus, while the flux loss in the hard band is still within 
the correlation scatter, the larger loss in the soft band throws the southern nucleus 
off the correlation of Eq. (|2.6|) . 



2.6 Effect of variability on the flux estimate 



IC342: variability of ULXs 

In the ASCA observations, the X-ray emission of IC 342 (a face-on spiral galaxy at 
3.9 Mpc) is powered by thre e main sources; two of them (source 1 and 2 according to 
Fabbiano fc Trinchierilll987l ) are ultraluminous X-ray binaries (ULX) while source 3 
is associated with the galactic cent re. Observations with higher angular resolution 
(ROSAT HRI. lBregman et al . 1993) showed that sources 1 and 2 are point-like while 
source 3 is resolved in at least three sources. Our main concern in determining the 
flux of this galaxy is the variability of the two ULX, which was assessed by a series 
of observations spanning seve ral years: IC342 was first ob served by Einstein in 1980 
f)Fabbiano fc Trinchieril Il987t) then by ROSA T in 1991 (iBregman et al.lll993h . by 



1998h and 2000 (|Kubota et alTl200ll) . and bv XMM 



ASCA in 1993 (jOkada et al 
Newton in 2001. 

In Table I2~T)1 we report soft X-ray fluxes for sources 1 and 2. We have chosen the 
soft X-ray band due to the limited energy band of both Einstein and ROSAT; the 
fluxes observed with these sa tellit es were obtained from the count rates reported in 
Fabbiano fc Trinchieri (1987) and Bregm an et al.l (|l993l) assuming the powerlaw and 
multicolor disk examined in Kubota et al.l ((2001 ) for source 1 and 2 respectively; we 



3 Assuming A v = 4.5 • 10 22 A^ H (Galactic value). 
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2 The X-ray luminosity as a SFR indicator 







Src 1 Flux 


Src 2 Flux 


Year 


Mission 


BB 


PO 


PO 


1980 


Einstein 


2.7 


4.1 


0.85 


1991 


ROSAT 


3.3 


2.7 


1.3 


1993 


ASCA 


16 




2.6 


2000 


ASCA 




5.0 


0.52 


2001 


XMM- Newton 




4.1 


1.3 



Table 2.3: 0.5-2.0 keV fluxes (in 10 13 erg s 1 cm 2 ) for sources 1 and 2 in 
IC 342 for a blackbody (bb) and a power-law (po) model. 



take ASCA 1993 and 2000 fluxes from lKubota et all (|200l|). XMM-Newton archival 
observations were reduced by us with SAS 5.3 and the latest calibrations available. 

Source 1 was in a low state (F . 5 _2 ~ 3 — 5 • 10~ 13 erg s _1 ) during the 1980, 1991, 
2000 and 2001 observations, and in a high state (^0.5-2 ~ 16 • 10~ 13 erg s _1 ) during 
the 1993 observation. The broad-band (0.5-10 keV) spectrum changed, its best-fit 
model being a disk black-body in 1993 and a power-law in 2000 and 2001. Source 2 
has also shown variability, its 0.5-2.0 keV flux oscillating between 0.52- 10~ 13 (ASCA 
2000) and 2.56- 10~ 13 (ASCA 1993) erg s -1 ; the main reason for this variability being 
the variations in the strongly absorbing column density, which was 9.9 ■ 10 21 cm -2 
in 1993 and 18 • 10 21 cm -2 in 2000. The spectrum was always a power-law. 

The high state for source 1 seems thus to be of short duration, and we feel confident 
that its time-averaged flux may be approximated with its low state flux. We thus 
choose to derive our flux estimate for IC342 from the ASCA 2000 observation, 
estimating the variation for the total flux of the galaxy caused by source 2 variability 
to be less than 10%. 



Variability in M82 



Hard (2-10 keV) X-r ay variability in M82 w as reported in two m onitoring campaigns 
with ASCA (in 1996. IPtak Griffith Jl 991 and RXTE (in 1 997. lRephaeli Gruber 
I20021 ). M82 was found in "high state" (i.e. 4 • 10" 11 < / 2 -io < 7 • 10" 11 ) in three 
out of nine observations with ASCA and in 4 out of 31 observations with RXTE. In 
all the other observations it was in a "low state" (2 • 10~ n < /2-10 ^ 3.5 • 10~ n ). 
A low flux level was also measured du ring the observations wit h other experiments : 



HEAO 1 in 1978 (|Orifiiths et al.lll979h : Einstein MPC in 1979 (IWatson et al. 
EXOSAT in 1983 and 1984; BBXRT in 1990; ASCA in 1993 dTsuru et al " 



1984); 



1997); 



BeppoSAX in 1997 (|Cappi et al.lll999i ): Chandra in 1999 and 2000, XMM-Newton in 
2001 . No variability was instead detected in the 0.5-2.0 keV band (IPtak fc Griffith^ 
1993). 

The high state has been of short duration: less than 50 days in 1996, when it was 
observed by ASCA, and less than four months in 1997. A monitoring campaign was 
also undertaken with Chandra, which observed M82 four times between September 
1999 and May 2000. We reduced the archival data, and found that the galaxy was 
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2.6 Effect of variability on the flux estimate 



always in a low state, with its flux slowly increasing from 1.2 • 10 11 to 2.9 • 10 11 



We do not attempt a detailed analysis of the variability (see iRephaeli fc Gruber 



2OO2]); however, we feel confident that, given the short duration of the high states 



and the fact that the difference between high- and low state flux is about a factor 2, 
the time-averaged flux of M82 can be approximated with its low state flux. We thus 
choose to derive our flux estimate for M82 from the BeppoSAX 1997 observation 
(/2-10 — 2.9 ■ 10 -11 erg s _1 cm -2 ), estimating the uncertainty caused by variability 
to be around 30%. 
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3 Star-forming galaxies in the 
Hubble Deep Field 



The X-ray and radio observations in deep fields reach limiting fluxes deep enough 
to detect star-forming galaxies at large redshifts z > 1. Thus they can be used to 
check whether the radio/X-ray relation holds also for distant galaxies. 

We consider the surveys performed in the Hubble Deep Field North (HDFN), a 
region of sky at a high Galactic latitude which was selected because of the very low 
extinction due to our Galaxy, centred at the (J2000) coordinates 12 h 36 m 49.5 s + 
62° 12' 58". The Chandra survey centred on the HDFN region (also called Chandra 
Deep Fie ld North, CDFN) w as initially performed for an assigned time of 1 million 
seconds (Brandt et al.N2001l ) and then extended with one other million secondes 



((Alexander et al.l l2003). It is customary to refer to the intermediate data release as 



the '1 Ms survey' and to the final one as the '2 Ms survey'. The final survey reaches 
a limiting flux of 2.5 ■ 10~ 17 erg s" 1 cm -2 in the 0.5-2.0 keV band for sources in the 
centre of the field. The HDFN wa s also imaged at radio wavelengths, with the Very 
Large Array (VLA) at 8.4 GHz (iRichards et alJll99«h and at 1.4 GHz ((Richards! 
a nd wi th the Westerbork Synthesis Radio Telescope (WSRT) at 1.4 GHz 
"l l2000h . The limiting fluxes for these radio surveys ar e all around 0.05 mJy . 




We searched for X-ray counterparts of radio sources in the lRichards et aL (1998) 
catalogue which contains optical and IR identifications allowing the selection of can- 
didate star-forming galaxies. Our selection criterium has been to include all galaxies 
with Spiral or Irre gular morphologies, k no wn redshifts and n o AGN signatures in 
their spectra (from Richards et al. 19981 or Cohen et al. l2000h. The X-ray counter- 
parts were initially searched for in the 1 Ms catalogue ((Brandt et al.ll2001[ ). When 
the final data release was made available, we reduced the X-ray data, generated 
a catalogue, and we looked for X-ray counterpars of the rad io sources. Our data 



reduction method followed rather closely the one described in lBrandt et al.l (|200lt ). 
Unless otherwise stated, all data reported here come from reduction of the complete 
2 Ms CDFN survey. 

The mean positional uncertainties of both Chandra (for on-axis sources) and the 
VLA are ~ 0.3", which added in quadrature give ~ 0.5". Using thi s value as the 



encirc ling radius for coordinate matching 5 galaxies were found in iBrandt et al. 
(2001). However, there are two effects that may increase this value: 



1. the shape and width of the Chandra PSF strongly depend on the off-axis and 
azimutal angles. Since the 2 Ms HDFN data consist of 20 observations with 
different pointing directions and position angles, there is no unique PSF model 
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3 Star-forming galaxies in the Hubble Deep Field 



Chandra 



VLA 



CXOHDFN 

CXOHDFN 

CXOHDFN 

CXOHDFN 

CXOHDFN 

CXOHDFN 

not present 

id. 

id. 

id. 



J123634.4+621212 (134) 
J123634.5+621241 (136) 
J123637.0+621134 (148) 
J123651. 1+621030 (188) 
J123708.3+621055 (246) 
J123716.3+621512 (278) 
in Brandt et al. (2001) 



3634+1212 
3634+1240 
3637+1135 
3651+1030 
3708+1056 
3716+1512 
3644+1249 
3701+1146 
3638+1116 
3652+1354 



Table 3.1: Identification of candidate star-for ming galaxies in thelRichards et al 
( 1998h catalogues with their X-ray counterpart ((Brandt et al.l l2001). Where appli- 
cable, the entry number in the Chandra catalogue is shown in parenthesis. 



even for near on-axis sources; 

2. a displacement between the brightest radio and X-ray positions, induced e.g. by 
an ultraluminous X-ray binary placed in a spiral arm and dominating the X-ray 
emission. 

Thus, by making cross-correlations with increasing encircling radii we found that 
the number of coincidences increases up to a radius of 1.0", yielding 7 matchings. 
There are no further coincidences up to a radius of several arcsecs, indicating that 
the sample should not be contamined by chance coincidences. 

By performing the same cross-correlation on the 2 Ms CDFN data, we added 
4 more galaxies with coincidences within 1.0", thus totalling 11 sources. One of 
these the 11 galaxies was dropped since it is probably an AGN on the basis of 
X-ray spectral and variability properties, as discussed below. The selected objects 
are listed in Ta ble 13.11 Fluxes at 1. 4 GHz and spec tral slopes were retrieved from 
Richardsl J2000) in 9 cases, and from [Garrett (2000) in one case. 



While the rest-frame 0.5-2.0 and 2.0-10 keV luminosities could be obtained by 
K-correcting the observer's frame counts, this would imply the assumption of a 
spectral shape; but none of the deep field galaxies is individually detected in the 
hard band, so that any constraint on their spectra obtained with the use of a hardness 
ratio diagram is too loose to be significant fFig. I3.1)) . However this problem may be 
partially circumvented by resizing the X-ray bands in the observer's frame according 
to the redshift of the objects - i.e. the counts extracted with this method correspond 
to the photons actually emitted in the rest-frame 0.5-2.0 and 2-8 keV bands. It is 
therefore possible to give better constraints to the spectra of the deep field galaxies 
and derive better estimates of their luminosities. 

Thus, we have redefined the soft and hard 1 bands as the [0.5; 2.0/(1 + z)] and 



1 Since Chandra has very poor sensitivity between 8 and 10 keV, the use of the reduced 2-8 keV 
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Figure 3 .1: Observer's-fram e hardness ratios for the 7 galaxies in the CDFN 
catalogue by lBrandt et al.l ((200 11 ): there are only upper limits since none is detected 
in the 2-8 keV band. The dotted lines show the loci for power-law spectra with 
slope T = 2.1 and different intrinsic rest-frame absorption; the dashed lines are for 
spectra of different slopes and no intrinsic absorption. 



[2.0/ (1 + z); 8/ (1 + z)] intervals, respectively. Another advantage of this procedure 
is that the higher the redshift, the more akin the new hard band is to the zone 
of maximum sensitivity of Chandra (~ 1 — 4 keV). Note that since the ACIS-I 
detector has almost no sensitivity below 0.5 keV, we fixed this energy as the lower 
limit for count extraction. For the two highest redshift galaxies, this reduces the 
soft band to 0.5-0.9 keV, still significantly larger than the ACIS-I energy resolution 
(FWHM < 100 eV). 

We extracted counts in circular regions (with radii of 5 pixel, i.e. ~ 2.5") around 
our selected targets; the background counts were taken in circular annuli surrounding 
the targets, with outer radii of 12 pixels (~ 6"), respectively. Counts and rates are 
reported in Table EOl Best-fit slopes reproducing the soft/hard count ratio were 
derived by assuming a power-law spectrum with Galactic absorption. We find that 
all objects have spectral slopes falling in the range V ~ 1.5-2.7 ( Table ET2*|) . To check 

band enhances the signal/noise ratio. Note that while counts are extracted in the 2-8 keV 
band, fluxes and luminosities are always extrapolated to the 2-10 keV band. 
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3 Star-forming galaxies in the Hubble Deep Field 
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Figure 3.2: The radio/X-ray luminosity relation for galaxies in the Hubble Deep 
Field (filled circles). Squares: local galaxies as in Fig. ([2.1)1 : dotted lines: linear fits 
for local galaxies (Eq. 12.6)1 . Besides the 10 deep field galaxies discussed in the text, 
here is shown also the point relative to the galaxy hd2-264.1 (in the Hubble Deep 



Field, at z = 0.475) which is the only galaxy from the sample of ICardiel et al. 



(200 



^) to have an X-ray detection (Lo.5-2 = 1.5 • 10 41 erg/s, L 2 _io = 5 • 10 40 erg/s, 



L 1A = 6 • 10 29 erg/s/Hz). 
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whether these spectra are consistent with those of the galaxies in the local sample 
we calculated the observed soft/hard flux ratio for galaxies in the local sample of 
Sect. 12.11 the median value for this flux ratio is 0.95 leading to a slope T ~ 2.1. The 
count ratios for each of the six deep field galaxies are consistent within l-2a with 
the T ~ 2.1 slope. 



One of the sources (#194 in iBrandt et alJ I2001L not shown in Table \A.2\ with 



z = 1.275) has an upper limit on the soft X-ray counts and a count ratio not 
consistent at the 5a level with the unabsorbed T ~ 2.1 spectrum: it requires an 
inverted spectrum (r < 0.1) if no absorption is assumed, otherwise, if we assume 
T = 2.1, the intrinsic absorbing column has to be iVn > 3.4 • 10 22 cm -2 . We also 
find that the flux of this source is strongly variable, since all detected counts were 
recorded in the first half of the survey, and no one in the second half which was 
performed about one year later. Thus, from its spectral and variability properties, 
we classified this source as an AGN and dropped it from the sample considered here. 

With the best-fit slopes we derived soft and hard band fluxes and luminosities 
(Fig. IH.2jl . The linear radio/X-ray correlations hold also for the Deep Field galaxies; 
the dispersion of the relations given in Sect J2. 21 is not changed by the inclusion of 
the deep field objects. 

A similar result was also found by Ba uer et al. (2002) by analyzing a sample of 
20 galaxies selected at 1.4 GHz in the Chandra Deep Field North. Their sample is 
partly overlapping with the one here described (not identical, since ours was selected 
at 8.4 GHz). 
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Fluxes and Luminosities: Deep Sample 



Soft X-rays Hard X-rays Radio 





Source 


z 


Net Cts. 


F 


L 


Net Cts. 


F 


L 


r 


F 


L 


a 


134 


3634+1212 


0.456 


55 


20 


13 


39 


28 


18 


1.9 


210 


1.4 


0.7 


136 


3634+1240 


1.219 


11 


5.6 


40 


44 


22 


160 


1.4 


180 


13 


0.7 


148 


3637+1135 


0.078 


27 


7.2 


0.10 




<4.6 


<0.065 




96 


0.014 


0.6 


188 


3651+1030 


0.410 


24 


12 


6.1 


4.3 


<2.8 


<1.4 




83 


0.42 


0.6 


246 


3708+1056 


0.423 


22 


7.5 


4.1 


13 


9.8 


5.3 


1.4 


36 


0.20 


0.4 


278 


3716+1512 


0.232 


22 


6.6 


0.94 


8.7 


8.2 


1.2 


1.4 


160 


0.23 


0.2 


new 


3644+1249 


0.557 


9.2 


3.4 


3.5 


7.3 


3.8 


3.9 


2.1 


31 


0.32 


0.7 


new 


3701+1146 


0.884 


6.7 


3.8 


12 


6.9 


2.8 


8.9 


3.4 


93 


3.0 


0.5 


new 


3638+1116 


1.018 


4.6 


2.8 


13 


7.4 


2.8 


13 


2.2 


92 


4.2 


1.1 


new 


3652+1354 


1.355 


3.8 


<2.4 


<23 


16 


6.2 


59 




20 


1.9 





Table 3.2: Data for deep field galaxies (counts, fluxes, luminosities). V is the best fit X - ray sl ope (photon index), a the 
radio slope (energy index). Sources are identifi ed via their entry num ber in the lBrandt et all ([20011 ) catalogue (cfr. Table l3~Tj) 
and their denomination in the VLA catalogue ((Richards et al. 1998| ). X-ray fluxes in 10 -17 erg s _1 cm -2 , radio fluxes in /xJy. 
X-ray luminosities in 10 40 erg s -1 , radio luminosities in 10 30 erg s -1 Hz -1 . X-ray counts are extracted in redshifted bands 
(soft band := 0.5-2.0/(1 + z) keV, hard band := 2.0/(1 + z)— 8/(1 + z) keV); the counts shown here are already background 
subtracted. X-ray fluxes and luminosities are in rest-frame 0.5-2.0 and 2.0-10 keV bands, radio ones at rest-frame 1.4 GHz. 
For source #194, the absorbed flux and luminosity are quoted. The unabsorbed luminosity is < 4 - 10 42 erg s _1 . Upper limits 
on fluxes and luminosities were assigned whenever the estimated error on the flux was of the same magnitude of the flux. 



4 The X-ray number counts and 
luminosity function of galaxies 



An estimate of the contribution of star-forming gal axies to the cosmic X-ray back- 

ground (XRB) has been a ttempted several times (e.g. Bookbinder et al.ll980l . Griffiths fc Padovani 
1990^ Moran et al. 1999I ). The main purpose for the earlier studies was the possi- 
bility to explain the flatness of the XRB spectrum via the X-ray binaries powering 
the X-ray emission of these galaxies (High Mass X-ray Binaries have flat spectra 
with slopes T ~ 1.2). Although AGN have since a long ti me been recognized to 
provide by far the m ost important contribution to the XRB ( Setti fc Woltier 1989; 



Comastri et al. 1995j ). the ongoing deep Chandra and XMM-Newton surveys offer 



unique opportunities to both test the AGN models and pin down the contribution 
from other kind of sources. 

A study of the X-ray number counts and of the luminosity function (LF) of galaxies 
is also valuable in the perspective of studying the evolution of the X-ray luminosity 
function (XLF) and the cosmic star formation history. In this chapter, the X-ray 
luminosity function and number counts of star forming galaxies are determined 
making use of different approaches, such as objects selection from X-ray surveys, 
conversion of radio and FIR LFs and number counts. An estimate is also performed 
of the contribution to the XRB. Finally, the possibilities to use X-ray observations 
to determine the cosmic star formation history are discussed. 



4.1 The X-ray luminosity function of galaxies 



The availability of multiw avelength photometry and spectral information f B arger et al 
2003; ISzokolv et al 2004 ) for the X-ray objects in the Chandra Deep Fields allows a 
direct identification of the star forming galaxies in these surveys. This is however a 
difficult task, for two reasons: i) at the limiting fluxes of current X-ray surveys, the 
majority of detected sources are AGN; ii) it may prove to be difficult, for objects at 
redshifts z ~ 1, to distinguish a bright starburst galaxy from a low luminosity type-2 
Seyfert galaxy (cfr. Sect. I1.5j) . Thus, in order to pick up normal galaxies in the X-ray 
survey samples, extreme care has to be put in the choice of the selection tools. We 
will not venture into a comp rehensive discussion of the s elect ion tools (good start - 
ing points for this mig ht be iMachalski fc Gondonl Il999l and iNorman et al.l 120041 ) , 
but just briefly illustrate a few of them that will be used throughout this chapter. 

A direct determination of the X-ray luminosi ty function (XLF) of galaxies has 
been recently attempted in INorman et al. (|2004h . where a sample was defined con- 
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4 The X-ray number counts and luminosity function of galaxies 



taining 210 galaxies w i th known redshift fro m the Chandra Deep Field catalogues 
( Alexander et all 120031 : iGiacconi etafll2002h . A Bayesian approach was chosen to 
derive a selection probability from the values of three different parameters: 

• the 0.5-2.0 keV X-ray luminosity for star forming galaxies is usually less than 
~ 10 42 erg/s; 

• star forming galaxies have a softer spectrum than AGN; this translates in 
selecting objects with an hardness ratio HR < —0.8 (with HR = (H—S) / (H+ 
S), where H and S represent the 2.0-10 and 0.5-2.0 keV fluxes respectively); 

• an X-ray/optical flux ratio Log(Fx_ ra y/-Popt) < — 1 (see below). 
The optical flux is defined as 

Log F opt = -OAR - 5.5 (4.1) 

were R is the red magnitude in the Kron-Cousins system. An X-ray/optical flux ratio 
equal (in logarithm) to —1 is usually taken as an a pproxim ate boundary betwee n 
normal galaxies and Seyferts ([Maccacaro et al.lll988t see also iMcHardv et!Hl2003h . 

We will consider the local differential luminosity function <£>(L)dLogL, which is the 
comoving number density of sources per logarithmic interval of lumi nosity. A binned 
lumin osity function (LF) was derived with the method develo ped inlPage &: Carrera 
(2000). which is a variant on the classical l/V m3iX method bv lSchmidtl (|l968l ). In a 
given bin of redshift and luminosity, the density of galaxies may be written as 



N 



(4.2) 



where dV/ dzdQ is the comoving volume per unit steradian sampled by the surveys. 
For a flat universe with non-zero cosmological constant we may take 



d^ 

d7 





' r dz' ' 


Go) E(z) 


[Jo E{z')\ 



(4.3) 



with 



E{z) = ^m(1 + ^) 3 + ^a- (4.4) 
The redshift distribution of the lNorman et al.l(j2004 ) sample is shown in Fig. (j4.1jl . 



Two redshift bins were considered, z < 0.5 and 0.5 < z < 1.2, with mean redshifts 
z = 0.27 and 0.79 respectively, in order to have a comparable number of galaxies in 
both bins. 

The XLF is shown in Fig. ([4.2)1 . iNorman et al.l (J2004J) also determined the evo- 



lution of the LF, by fitting the XLFs in the two redshift bins with two power laws. 
Considering a pure luminosity evolution, they found L(z)/L(0) = (1 + z) 2 ' 7 . This 
result will be checked in the following sections against FIR- and radio-based deter- 
minations of the XLF. 
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4.2 Determining the XLF from the FIR and radio luminosity functions 
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Figure 4.1: Redshift distribution for 
the Bayesian sample of X-ray selected 



norma l galaxies defined in lNorman et al 



(|2004f ). Dotted and dashed line: redshift 
distributions for galaxies in the Chan- 
dra Deep Fields South and North, respec- 
tively. Continous line: sum of the North 
and South subsamples. 



Figure 4.2: X-ray luminosity func- 
tion fo r the galaxies in the B ayesian sam- 
ple bv iNorman et al.l ((2004J) . Triangles: 
LF of objects with z < 0.5; circles: LF 
of objcts with z > 0.5. The median red- 
shift of the two bins are z = 0.27 and 0.79 
respectively, thus the two data sets may 
be regarded as the XLFs at the median 
redshifts; the evolution of the LF from 
z = 0.27 to z = 0.79 is clearly visible. 



4.2 Determining the XLF from the FIR and radio 
luminosity functions 

Infrared surveys are a powerful method to select star forming, spiral galaxies, since 
the bulk of the far and near infrared emission is due to reproc essed light from star 



Franceschini et al. 2001 



formation, with AGN r e presenting only a minor population (|de Jong et al 

Elb az et al 



1984; 



2002). In the following, we summarize the 



current determinations of the FIR luminosity functions, which are then converted 
to the X-rays. 

We will only consider a pure density evolution of the form (1 + z) v , which is the 
simplest and most used form by the authors of the papers from which we draw the 
FIR LFs. This is equivalent, in a more rigorous formalism, to a bivariate LF of the 
form (fb(L,z) dLogLdz with the assumption that the LF is separable for L and z, 
so that we may w r ite <Pb (L, z) = ip(L)(l + z) v . 

Sa unders et al. (1990) defined a sample of 2,818 galaxies matching common se- 
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4 The X-ray number counts and luminosity function of galaxies 
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Figure 4. 3: Left panel: Redsh ift distribution for galaxies in the IRAS PSCz 
survey, from Saunders et al. ( 2000l ) . Right panel: Luminosity-redshift relation for 
galaxies in the ELAIS survey, from lSerjeant et al. (2004). 



lection criteria from different IRAS samples, with a flux limit around 0.6 mJy at 
60/x, completeness at 98% level, redshifts z < 0.05, and derived a 60/i luminosity 
function. The best-fit model has the shape 



L 



l-Q 



2(T 2 



Log? 



<P(L) =tf \ — ) <>xp 

with parameters at present epoch 

(p* = (2.6 ± 0.8) ■ 10~ 2 h 3 Mpc" 3 
a = 0.72 ± 0.03 and L 



Mpc" 



(4.5) 



a = 1.09 ±0.12, 

+3.06\ 



(4.6) 



(2.951?;°?) -10 8 h- 2 L Q 

(shown as the dotted curve in Fig. 14. 4 j) and the evolution (parameterized as pure 
density evolution) 

(l + zf 7±2 -\ (4.7) 

The 60/x luminosity function was revised by Ta keuchi et al. (2003) by enlarging 
the galaxy sample. 15,411 galaxies from the Point Source Catalog Redshift (PSCz) 
were used, covering 84% of the sky with a flux limit of 0.6 mJy at 60/i. Using the 
same parameterization of Eq. (|Oj) r iTaWhi et all ([iool found 



(p* = (2.60 ± 0.30) ■ 10~ 2 h 6 Mpc" d , a = 1.23 ± 0.04, 
a = 0.724 ± 0.01 and L* = (4.4 ± 0.9) • 10 8 h~ 2 L Q . 



In Fig. 14.41 it is shown as th e dotted curve. They also found that the (1 + z) 6 7 
density evolution proposed by ISaunders et al. (|l990i) is not consistent with their 



lr T. Takeuchi (priv. comm.) recently found an error in the value of L* reported in his paper. The 
one reported here is the correct one. 
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Figure 4.4: IRAS and ISO local 
(i.e. z = 0) luminosity f unctions. Solid 



curve: IRAS 60fi LF from lTakeuchi et al 
(l200aj). Dotted cur v e: IR AS 60// LF 



Saunders et al.l (jl990h. The two 



from 

IRAS LFs coincide at bright luminosi- 
ties. Dashed c urve: ISO 90/i LF from 
Serjeant et al.l (J2004J); the shaded area 
corresponds to the la errorbar. While 
the IRAS LFs are well constrained in the 
whole range of luminosities shown here, 
the ISO LF suffers from incompleteness 
for L < 1O 9 L . The monocromatic lu- 
minosities have been converted to bolo- 
metric ones by assuming vL v = cost. 
H = 70 is assumed. 



Figure 4.5: IRAS, ISO, and ra- 
dio luminosity functions converted to 
the X-r ays. Solid curye: I RAS 60/i 



LF from iTakeuchi et "all ll2003h. Dotted 



curve: IRAS 60/i LF from lSaunders et al 
1990 h- Poin ts: ISO 90/i LF from 



Serjeant et al.l fl2004h . Dashed curve: 



radio LF 
( 2000h . 



from iMachalski fc Godlowski 



sample. A milder evolution (1 + z) 3A is reported to be the best-fit description of the 
data. Note however that the redshift distribution of the PSCz galaxies only span a 
very limited range (Fig. 14.3)) . thus the evolution exponent is very poorly constrained; 
this may also explain the factor-of-two difference in evolution between Takeuchi's 
and Saunders' results. 



Th e LF a t 90/i was recently determined by the ELAIS survey team ([Serjeant et al 



2001L |2004) with the PHOT and CAM i nstruments onboard the ISO satellite. It 



has not been fitted to a parametric form. iSerjeant et ajj ((20041) also found that the 
evolution, if parameterized as a (1 +z) v pure density form for z < 1 and no evolution 
beyond z ~ 1, has an exponent 2.4 < i] < 4.6 at 68% confidence; in the following we 
will assume rj = 3. The ELAIS LF is shown in Fig. 14.41 as the dashed curve, along 
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4 The X-ray number counts and luminosity function of galaxies 
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Figure 4.6: Comparison of the X-ray luminosity function with the FIR- and 
radio-derived ones. Symbols as in Figs. f)4.2l 14.5)) . The FIR- and radio-based LFs 
have been evolved (see text) to the median redshifts of the two bins in which the XLF 
was calculated (i.e. z = 0.27 in the left panel, z = 0.79 in the right panel). For the 
Saunders et al.l (|l990i ) IRAS LF, the milder evolution determined bv lTakeuchi et aL 
(200^) was used. 
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4.2 Determining the XLF from the FIR and radio luminosity functions 



with la errors. 



Conversion of the LF to the X-rays 



The FIR luminos ity functions may be conve rted to th e X-rays using the approach 
first developed in Avni fc Tananbaum ( 1986[ ) (see also Georgantopoulos et al. 19991 
and lNorman et al.ll2004l ). 



-t-oo 

V9 X (Log L x ) = J v9 60A1 (Log L 60ai ) P (Log L x |Log L 60m ) dLog L 60m 



(4.9) 



where P (Log Lx|Log ^60^) is the probability distribution for observing Lx for a 
given L 60m . Considering the radio/FIR relationship of Sect 12 .2\ this may be taken 
as a gaussian, yielding (for the 0.5-2 keV band) 



P(Log L . 5 -2|Log L 60/Lt ) 



L °g - L 60M+ 9 ' 05 ~ Log L 0-5-2 



2na 



(4.10) 



where a ~ 0.30 as estimated in Sect. 12.21 To convert the ELAIS LF, a Sqq^/ = 
0.66 flux ratio may be derived by assuming vL v = cost. The same figure can be 
obtained by considering a grey body law (Eq. II. 8|) with f3 = 1 and temperature 
T = 40 K. The converted luminosity functions are shown in Fig. 14.51 



The XLF from the radio LF 

Information from the radio surveys may also be used to derive an XLF. We con- 
sider _^]£_^airrple_of_j^dio detected galaxies in the Las Campanas Redshift Sur- 
vey ( Machalski fc Condon 1999h . which comprises 1157 galaxies with magnitude 
R < 18.0 and detected in the NRAO VLA Sky Survey (NVSS), i.e. brighter than 
2.5 mJy beam" 1 at 1.4 GHz. These galaxies were classified as "starburst" (502 
galaxies) or "AGN" (655 galaxies) on the basis of t heir FIR-radio flux ratio s, 25m- 
60/i spectral indices and radio-optical flux ratios. Machalski k, Godlowski (j2000) 
derived the LFs for both classes of objects; here we consider only the "starbursts" . 

For a LF with the functional form of Eq. (J4.5|) . Machalski fc Godlowski (2000) 
give the set of parameters 



= (7.9; 



-3.8\ 



■ 10" d /^ Mpc" d dex" 



a 



1.22 ±0.27, 



a = 0.61 ± 0.05 and Log L* = (28.13 ± 0.37) h~ 2 erg s" 1 Hz" 1 . (4.11) 

We converted the radio LF into an XLF following the same method as in the 
previous paragraph, and considering the radio/X-ray relation of Sect. 12.21 The 
resulting XLF is shown in Fig. (J4.5|) along with the ones derived from the infrared 
surveys. The radio-based LF agrees well with the FIR-based ones for Lx ^ 10 40 
erg/s, while it lies a factor of ~ 2 below the FIR-based LFs at lower luminosities. 
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4 The X-ray number counts and luminosity function of galaxies 



The redshifts of the Las Campanas "starburst" galaxies are in the range 0.01 < 
z < 0.075. As in the case of the IRAS galaxies, this redshift range is too small to 
all ow a reliable estimation of the e volution. The constraints on the evolution found 
by Machalski fc Godlowski (jio OO) are too loose to either confirm of reject cosmic 
evolution for the Las Campanas galaxies. 



4.3 Comparison with the XLF 



To compare the pure luminosity evolution determined bv lNorman et al.l 1)20041 ) with 
the pure density forms used above, we assume that the XLF may be well represented 
by a power law whose slope remains constant with varying redshifts. We remind, 
in fact, that for a power law shaped LF the two forms of evolution cannot be 
distinguished. In fact, it may easily be shown that a LF with pure luminosity 
evolution with the shape 



(p L(z) 

transforms into a pure density form as 



[L /(i+ z y] 



^ = ^(1 + 2)%' 



(4.12) 



(4.13) 



Since the XLF derived bv iNorman et aL has a ~ 1.1, their (1 + z) 2 ' 1 lumi- 

nosity evolution is equivalent to a (1 + z) 2 ' 9 density evolution — which agrees well 
with the evolution exponents found for the FIR LF s. 

A comparison between the Norman et al.l ( 2004 ) LF and the FIR/radio derived 
ones is made in Fig. ()4.6j) . To make the comparison the easiest possible, the 
FIR/radio derived XLFs were evolved with th eir own evoluti o n exp onents, up to the 
mean redshifts of the two bins in which the iNorman et al.l (|2004h was calculated. 
We found that the X-ray derived LF is in agreement with the FIR-radio derived 



41.5> 



LFs for Lx < 10 , but it is significantly flatter at its bright end (Lx > 10 
where it has a density of galaxies which is one order of magnitude larger than that 
of the FIR-radio galaxies. This may b e symptomatic of a r esidual fraction of low 



luminosity AGN in the sample used bv lNorman et al . (2004). 



4.4 The observed X-ray Log TV-Log S of galaxies 



We turn now to an analysis of the X-ray number counts of galaxies, since they 
represent the most direct testbed for both the selection criteria for normal galaxies 
in X-ray surveys, and for the shape and evolution of different LFs. 



Firs t we considered the Log iV-Log 5* of the sample defined in INorman et al. 
(2004J) and discussed in Sect. 14.11 which is shown in Fig. (j4.7J) as the long-dashed 
histogram. 

Then we considered a sample defined in IBauer et al. (|2002h . who investigated 
the link between the faint X-ray and radio sources in the Chandra Deep Field North 
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Figure 4.7: X-ray counts in the Chandra Deep Fields. The thick, upper line 
and the horn-shape d symbol show the observed Log iV-Log S for all X-ray sources 
i n the Deep Fields ( Moretti et al. 2003| ) and the limits from the fluctuation analysis 
( Mivaji fc Griffiths! 2002a Ibfh" The histograms show the Log iV-Log S for normal 
galaxies. Short-dashed (lower) histogram: sources with Log (X-ray/optical flux 
ratio) < — 1. Dotted (middle) histogram: X-ray sources with a radio detection an d 
an optical spectrum compatible with ion excition by H II regions dBauer et al . 2002) 



Long- dashed (upper) histogram: sources from the Bayesian sample of iNorman et al 
f)2004h . 
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(CDFN); 11 sources from their catalogue could be reliably selected as star formation 
powered systems, on the basis of the reported optical spectrum (i.e. classified as 
'emission l ine galaxies" ) . By cons idering the sky coverage of the CDFN survey, as 



reported in Alexander et al. ( 200.^ 1. we derived the Log iV-Log S for these sources, 
which is shown in Fig. (|4.7jl as the dotted hist ogram. 



We also considered 97 CDFN sources from Barger et al. ( 2003| ) . which were se- 



lected for having an X-ray/optical flux ratio < — 1. Although this selection criterium 
may be regarded as somewhat crude, it is indeed useful since it allows an immedi- 
ate object selection with the least possible number of parameters. The resulting 
Log A^-Log S is plotted as the short-dashed histogram in Fig. ()4.7|) . 

We found that there is a quite good agreement (within a factor of 2-3) between 
the three observed Log A^-Log S described above. In the following sections, the 
X-ray Log A^-Log S will be checked against predictions made on the basis of the 
radio counts for the sub-mJy population and from the FIR and radio LFs. 



4.5 X-ray number counts from radio surveys 



We consider the radio sub-mJy populat ion associated with fain t blue galaxies at high 
redshifts (22 < V < 27, 0.5 < z < 1.5 : IWindhorst et al.lll990l) representing an early 
era of star f ormation in the universe ((Haarsma et al.ll2000f ). This strongly evolving 
population (IFranceschini et aT1l200lh accou nts for the majority of the number counts 



below ~ 0.5 mJy ((Windhorst et al 



- 19851) and contrib utes about half of the radio 

cosmic background at 1.4 CHz (lliaarsma et al. 2000[ ). The deepest radio sur veys 
have bee n performed at 1.4 GH z ((Richards 20Q0J), 5 GHz ( Fomalont et al. 1991 ) and 



8.4 GHz ( Richards et al. 1998| ). In order to derive the X-ray counts from the radio 
ones a full knowledge of the redshift distribution and spectra of the sources would 
be required. Under the simplifying assumption that the sub-mJy population lies at 
a redshift z, so that the K-correction term is the same for all sources, the differential 
counts are obtained as 

n(S x ) = n(S 1A ■ «* 4 ) (4.14) 

where n(S) are the differential number counts, and ny~ A = ky a (z, a,T) is the X-ray 
band/1.4 GHz luminosity ratio (eqs l2.6p2.10j) which depends on the redshift and on 
the radio and X-ray spectral indices via the K-correction: 



nf A (z,a,T) 



K 



X I 
1.4 



1 + _ r( r-i )+Q 



(4.15) 



where, as customary, a denotes the radio spectral index ('energy index', defined by 
S(u) oc v~ a ) and Y is the X-ray spectral in dex ('phot o n ind ex 7 , S{u) oc z/ 1_r ). 

Concerning the sub-mJy Log A^-Log S, iRichards! (|2000i) gives n{S) = (2.51 ± 
0.13) ■ io _3 5' _2 ' 38±0 ' 13 deg~ 2 Jy _1 as a best- fit to the differential n umber counts 
at 1.4 GHz in the range 45-1000/iJy, while iFomalont et all (|l99l[ ) find n(S) = 
1.2 ■ io- 3 S , - 2 - la±0 - 19 at 5 GHZ in the range 16-1500 /xJy. The number density at 4 
/j,Jy (as estimated from fluctuation analysis) is consistent with extrapolation of the 
16-1500 fiJy slope. 
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The mean radio spectral index is in the range a = 0.3-0.7. iFomalont et al.l (J19_9jJ) 
report their distribution as having a moda of a = 0.5, a median of a = 0.38, 
and an average of a = 0.28 (indices measured between 1.5 and 5 GHz). From 
Richardsl ( 2000[ ) data we find an average a = 0.47 (between 1.4 and 8.4 GHz), 
when we consider detections at both frequencies, and a = 0.67 when we consider 
only detections at 1.4 GHz and treat upper-limits at 8.4 GHz as detections. The 
latter, steeper slope is also more consistent with the average index of our HDFN 
sample (Table l3~2|) . Here we assume a = 0.5, and estimate the uncertainty in the 
K-correction due to the radio spectral index to be around a 15 — 20%. 



According to the results of Sect 01 we further assume an average X-ray spectral 
slope T = 2.1. For the sake of simplicity we also assume that the objects are 
placed at z ~ 1, the m ean redshift for the sub-mJy galaxies which, according to 
Windhorst et al. ( 1990[ ). are distributed in the redshift interval 0.5-1.5 with a peak 



at z ~ 1. To estimate the effect on the counts due to the actual distribution of 
the sources in this redshift interval we consider a simplified case in which they are 
equally distributed at redshifts z = 0.5, 1.0 and 1.5. Since the effect enters in the 
computations via the K-correction term, we find that the predicted number counts 
would increase by only ~ 10%. As a last remark we notice that a fraction, as yet 
undefined, of the X-ray spectra might steepen at energies > 10 keV, thus entailing 
a decrease in the predicted source counts; however, given the redshift range under 
consideration, the predicted soft X-ray source counts should not be affected. 



The X-ray number c o unts predi cted with the above a ssumptions (i.e. a = 0.5, 
T = 2.1) from iRichardsj JlOOO) and Fomalont et al. I (Il99lh Lo g A— Log S are shown 
in Fig. ()4.8|) . along w ith observed counts fromlMoretti et al. ( 20031 ) and limits from 
fluctuation analysis (|Mivaii fc Griffiths! l2002al ]bh. The extrapolations of the radio 
Log iV-Log S below ~ 50/^Jy and down to a few /iJy at 1.4 GHz (~ 6 ■ 10~ 18 erg 
s" 1 cm" 2 in the 0.5-2.0 keV band ) are consistent with the limits from fluctuation 
analysis ( Mivaji fc Griffiths! l2002b| ) and do not exceed the X-ray expected number 
counts. It should be noted, however, that eqs. ()2.6|2.10|) may not ap ply to the entire 



sub-m Jy population. For instance, a fraction (up to a 15 — 20%, lHaarsma et al. 



200o]) of the sources may still present a sizable contribution from an AGN. It follows 



that the derived Log iV-Log S should be regarded as an upper limit to the X-ray 
counts from star-forming galaxies. 



As a final remark it should be pointed out that our results depend on the basic 
assumption of a strict linearity between the radio and X-ray luminosities. Had we 
assumed a non-linear relationship, such as Eq. (J2.4|) for the soft X-rays or Eq. (|2.8|) 
for the 2-10 keV band, we would have found an increase or decrease, respectively, of 
about 50% in the predicted counts at a flux level around 10~ 17 erg s" 1 cm" 2 . 
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Figure 4.8: X-ray counts derived from deep radio Log iV-Log S in the 0.5-2.0 (upper 
panel) and 2.0-10 keV (lower panel) bands. Histogra ms as in previou s figures. The short- 
dashed lines represent the 1.4 GHz Log iV-Log S (Ri chards! 12000) with K-corrections 
as suming z = 1, T = 2.1 and a = 0.5. The long-dashed lines represent the Log iV-Log S 
bv lFomalont et al.l (j.1991) reported at 1.4 GHz with a = 0.5 (K-corrections like for the 1.4 
GHz one). Thick line and horn-shaped symbols: total X-ray number counts and resuls 
from fluctuation analysis, respectively (see Fig. 14.7ft . The vertical dotted lines show the 
limiting sensitivities for the radio surveys. For easier reading, the horizontal scale is shown 
both in radio and X-ray fluxes, assuming the conversion factors of Sect. 12.21 
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4.6 X-ray counts from integration of the FIR and radio LFs 




Fig ure 4.9 : X-r ay counts derived from the IRAS luminosity function by 
Saunders et al. (|l990j ). Left panel: (1 + z) 6 - 7 density evolution; right panel: (l + z) 3A 
density evolution (see text). Thick line, horn shaped symbol and histograms as in 
the previous figures. The three thin curves, which converge at fluxes > 10~ 16 erg 
s _1 cm -2 , show the number counts from the integral of the LF; lower curve: the LF 
was integrated with z max = 1.1 and evolution as described in the text; upper curve: 



2; middle curve: z n 



2 but evolution stopped at z 
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Figure 4.10: X-ray counts de- 
rived from the I RAS l uminosity function 



Figure 4.11: X-ray counts derived 
from the ELA I S lum inosity function by 



Tak euchi et al. Symbols as in ISerieant et al.l Symbols 
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Fig. (gU). 



Fig. (|OJ). 
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4.6 X-ray counts from integration of the FIR and 
radio LFs 

The observed X-ray counts may also be checked against the number counts obtained 
by integrating the FIR- and radio based luminosity functions. However, it should be 
noted that the galaxies considered in the IRAS and ELAIS surveys lie at redshifts 
much lower than the Chandra Deep Fields ones. The redshift distributions for the 
PSCz and ELAIS surveys (Fig. 14.3)1 show that the bulk of the galaxies considered 
in these surveys is placed at z < 0.05 (PSCz) or z < 0.3 (ELAIS), while the X-ray 
faint galaxies are usually placed at z ~ 1 (cfr. Chap. El and Sect. 14.5)1 . Thus a 
comparison of the X-ray counts with the FIR-derived ones involves an extrapolation 
of the luminosity function at higher redshifts. Although this introduces a further 
source of uncertainty, it may be regarded as a possibility to test for the maximum 
sustainable evolution. 

The number counts are given by: 

2max ^max 

N{>S)= J J (p{L,z)^- dLdz (4.16) 

Zmin L m i n (2) 

where dV/ dz is the comoving volume between z and z + dz. In the following the 
integration will performed in the luminosity interval 10 39 < Lx < 10 42 erg/s and 
for two different intervals in 2: < z < 1.1 (cfr. Sect. 14.1)1 and over the enlarged 
redshift range < z < 2. We also considered as an intermediate possibility to stop 
the evolution at z = 1. 

The Log iV-Log S resulting from the integration of the infrared luminosity func- 
tions are shown in Fig.s (J4~TU|) . and (j4~TTJ). 



In t he left panel of Fig ()4.9j) we plot the counts obtained from the ISaunders et al 



(1990) ^F with the evolution (1 + z) 6 7 . It is immediately clear that this evolution is 



too strong, since the predicted counts lie much above the observed Log iV-Log 5" at 
flux es < 4 ■ 10~ 17 erg s" 1 cm -2 . If we assume a milder evolution (e.g. the one found 
by Takeuchi et al. 20031 (1 + z) 3 ' 4 ), then the predicted counts are consistent with 



the observed ones. In this case, a good agreement is also found with the predictions 
from the radio Log iV-Log S of Sect. 14.51 (cfr. Fig. 14.81 upper panel). 

Integrating the EL AIS LF gives number counts which are very similar to the 
Takeuchi et aD j2003) case (Fig. 14.11)1 The X-ray number counts obtained by in- 



tegrating the Las Ca mpanas LF (not shown) are also very similar to the IRAS 
( Takeuchi et~aL1l2003h and ELAIS LFs (for the radio LF we considered a pure den- 



sity evolution with (1 + z) 3 , i.e. the same evolution of the ELAIS galaxies). 



4.7 Final remarks on the X-ray counts 

We have derived estimates for the X-ray number counts at faint fluxes (5 ■ 10~ 18 < 
Sx < 10~ 15 erg s^ 1 cm~ 2 ) with various methods (selection of galaxies in X-ray deep 
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surveys, conversion of FIR and radio LFs, conversion of radio counts). We note that 
all of them agree in two points: 

• the fraction of normal galaxies at the current flux limit (~ 5 ■ 10~ 17 erg s _1 
cm _2 in the 0.5-2.0 keV band) of the Chandra Deep Field surveys is about the 
20% of the total number of objects; 

the X-ray number counts of normal galaxies shoul d overcome the counts from 



AGN at fluxes below 1-2-10" 17 erg s" 1 cm" 2 . Since iMivaii fc Griffiths! J2002b; 



suggested the emergence of a new population (beyond that of AGN which 
dominates at brighter fluxes) in the 0.5-2.0 keV band at fluxes around 10~ 17 
erg s _1 cm" 2 , it is tempting to identify this new population with the sub-mJy 
galaxies. 

We also note that it is unlikely that the X-ray Log iV-Log S of normal galaxies 
could sustain its slope much below 2-3- 10 -18 erg s _1 cm -2 , for two reasons: i) be- 
cause of radio /X-r ay correlation, the in tegrated radio emission from weak sources 
would diverge (cfr. iFomalont et al. 1991 ): ii) a much stronger evolution of the lumi- 
nosity functions would be required for the counts not to flatten below ~ 10~ 18 erg 
s -1 cm" 2 ; but this would be at odds with the limits from the fluctuation analysis 



( Mivaii fc Griffith jl2002bh . 



4.8 The contribution to the X-ray background 

The derived X-ray number counts can be integrated to estimate the contribution to 
the XRB. For the observe d 2-10 keV backgr ound we take the XMM-Newton value of 
2.15 • 10~ n erg s" 1 cm " " 2 (|Lumb et al.ll2002h. which is comp rised between the ASCA 
( Gendreau et al. 19951 ) and BeppoSAX ( Vecchi et aL 19991) figures . The integration 
of counts derived from the 1.4 GHz Log iV-Log S ( Richard's] 2000[ ) . performed in its 



validity range (5.9 • 10~ 17 — 1.3 • 10~ 15 erg s _1 cm -2 , corresponding to 45-1000 /iJy 
at 1.4 GHz), yields a contribution to the XRB of 4.0 • 
of the observed background). By extrapolating to 10~ 18 



10 13 erg cm 2 s 1 



contribution would increase to 2.3 • 10 12 erg cm 2 s ] 



(< 2% 

erg s _1 cm 2 (~ 1/xJy) the 
11%). Integration of counts 

from the flatter Fomalont et al. ( 199ll ) Log iV-Log S in the 1-1000 /iJy range yields 
a contribution of 1.4- 10~ 12 erg cm -2 s" 1 (6.4%). We also note that 1 /iJy at 1.4 GHz 
is a limit ~ 3 — 5 times fainter than the constraint from radio fluctuation analysis, 
and it is unlikely that the radio Log iV-Log S could sustain its slope below this limit 
(cfr. previous section). 



4.9 An X-ray estimate of the cosmic star formation 
history 

The current determinations of the cosmic star formation rate density (CS FRD) show 
a rapid increase of the size of one order of magnitude from z = to z ~ 1 (|Lillv et al 
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Figure 4.12: Esti mates of the cosmic star formation history based on different 



The black squares represent the X-ray 

0.2 7 and 
( 1995h 




estimates derived in iNorman et al 

0.79. The triangles r epresent the Ha SFR values from iGallego et al 



at the median redshifts z 

r 



at z ~ 0, iGronwalll (|1999h at z~ 0.05, and 
empty square represents the UV-selected Ha 



HoDkins et al. (2000) at z~ 1.3. The 
Pascual et al.1 (l200ll) value at z~ 0.24. 

| 2002h . The star gives 



The filled circles show the SFR densities f rom iTresse et all .. _ 
the UV-selected z=0.15 SFR density from Sullivan et al.l (|2000l ). The open circles 
give 2800A CFRS points from Lilly et al. ( 1996f ). without dust correction. We also 
plot the SFR history based on the 60 //m LF luminosity density at z ~ 0, including 



(T+z) 3 evolution. Th e 60 um luminos i ty den sity was computed from lSaunders et al 
( 1990L solid line) and lTakeuchi et al.1 (|2003l dotted line). 



1996J; Madau et al. 1996| ). At larger redshifts, either a peak at z ~ 1-2 and a slow 



decrease, or a plateau are found, if UV and Ha data are taken at their face value, 
or if they are corrected for extinction, respectively 

The main advantage in estimating the CSFRD from X-ray surveys would be that 
the X-ray emission from star forming galaxies does not suffer from absorption, so 
that one of the main sources of uncertainty in the CSFRD determinations is removed. 
However, the current X-ray surveys do not detect galaxies beyond redshifts ~ 1.2, 
nor any deeper survey is currently planned. Thus the most interesting part of 
the cosmic time (z ~ 3-5) is precluded to these studies. Nonetheless, it is still a 
good exercise to check if an X-ray estimate of the CSFRD is consistent with the 
determinations from data taken at other wavelengths. 

A determination of the cosmic star formation rate in a redshift interval can be 
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obtained by (cfr. eq. 12.11}) 



SFR(^) = 2.2 ■ 10 



-40 



<p{L)L dL. 



(4.17) 



In Fig. (14.12 ] ) we show the resu lts for the cosmic SFR at the mean redshift consid- 
ered for thelNorman et al.l (J2004J) sample, along with compilation of SFR estimates 
from Tresse et al.l (|2002j ). The overall agreement of the X-ray derived points with 
the ones estimated from infrared and Ha observations is very good. This may not 
be surprising, because the X-ray and FIR LFs agree quite well and the CSFRD con- 
veys the same information contained in a luminosity function, but this result indeed 
confirms that X-ray observations may successfully be used to derive an estimate of 
the CSFRD and of its evolution. 



65 



4 The X-ray number counts and luminosity function of galaxies 



66 



5 Metallicity enhancement in 
starbursts: the case of M82 



The signature of the star formation (SF) history of a galaxy is imprinted in the 
abundance patterns of its stars and gas. Determining the abundance of key elements 
released in the interstellar medium (ISM) by stars with different mass progenitors 
and hence on different time scales, will thus have a strong as trophysical impac t in 
drawing the global picture of galaxy formation and evolution (|Mc William! 1 1 9 9 7) . 

Determining the abundances of SB galaxies is not only relevant to constrain their 
SF histo ry but also offers the u nique chance of directly witnessing the enrichment of 
the ISM flMaeder fc Contill994h . Metals locked into stars give a picture of the enrich- 
ment just prior to the last burst of SF, while the hot gas heated by SNe II explosions 
and emitting in the X-rays should trace the enrichment by the new generation of 
stars. Nebular metallicities from the cold gas should potentially trace the stellar ones 
if the c ooling and mixing ti mescales are slow (up to 1 Gyr), as suggested by some 
models ( Tenorio-Taglel 1996| ) since the gas enriched by the ongoing burst is still too 
hot to be detected in the optical lines. On the contr ary, if rapid (a few Myr) cool- 
ing a nd mixing occur as suggested by other models ( Recchi. Matteucci fc D'Ercold 
2001), the nebular abundances could trace the gas enriched by the new population of 



stars (similarly to X-rays). Hence, important constraints on the cooling and mixing 
timescales of the gas can be provided by comparing the metallicity inferred from 
stars and hot plasma with those of the cold gas. 

So far, the metallicities of distant starburst (SB) galaxies have be en mainly derived 



by measuring the nebular lines associated to their giant HII regions (ptorchi-Bergmann. Calzetti fc Kinnev 



1994 ; Coziol et al. 19991 : Henry. Edmunds &: Koppenl 2000l ) . Direct determination of 



O and N abundance s can be obtained once electron temperature (T e ) is well estab- 
lished (jMcCalllll984l ). At high metallicities T e decreases due to strong lines cooling 
and its major diagnostic (the weak [OIII] A4636) is no longer observable. Several 
alternative methods based on other strong lines were pro posed in the past years 
fjPagel et alJll97d : iDooita fc EvanilllQSd : iMcGauerhl Il99lh to derive T e , but all of 
them rely on a large number of physical parameters a nd model assump tions for the 
geometry and the nature of the ionizing populations (|Stasinskall2001[ ). Depending 
on the adopted calibrations, quite different abundance sets can reproduce the ob- 
served line ratios. Moreover, stellar abundances can be poorly constrained from 
optical spectra, since the nebular emission strongly dilutes the absorption lines and 
dust can heavily obscure the central regions where most of the burst activity is 
concentrated. 

However, abundances in SB galaxies can be also measured from absorption stellar 
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features in the near IR and/or emission lines in the hot X-ray gas. In the near IR the 
stellar continuum due to red su pergiants (RSGs) usua l ly largely dominates o ver the 
possible gas and dust emission (|01iva fc Origlialll998t lOridia fc Olivall2000h . while 
the X-ray emission from SB galaxies at energies lower than ~ 2 keV is mainly due 



to hot plasma heated by SN explosions (jDahlem et al.lll9 98) (Sect. fT3|) 



In order to test the above discussed enrichment scenario we have started a pilot 
project to measure the metallicity enhancement in a sample of starburst galaxies, 
for which we obtained high resolution infrared (J and H band) spectra with the 3.6 
m Italian Telescopio Nazionale Galileo (TNG) and archival data from the XMM- 
Newton and Chandra missions. Our sample comprises M82, NGC253, NGC4449 and 
the Antennae, sampling a difference of two orders of magnitude in star formation, 
as it ranges from the 0.3 M /yr of NGC4449 to the 30 M Q /yr of the Antennae. 

The preliminary results which have been achieved for M82 with the available 
XMM-Newton archival data (a 20 ks exposure which provides high resolution data, 
although with a low signal to noise ratio) will be presented throughout this chapter. 
Sect. briefly reviews the nebular abundance estimates of M82, as published so far 
in the literature. Sect. 15 .21 shows our near IR spectra of the nuclear region of M82 and 
the derived stellar abundances. Sect. I5.3l presents a re-analysis of the XMM-Newton 
and Chandra nuclear X-ray spectra and the hot gas-phase abundances. Sect. 15.41 
summarized the overall chemical enrichment scenario in M82 as traced by the stellar, 
cold and hot gas-phase abundances. 



5.1 Nebular abundances of M82 from the literature 



M82 is a prototype of SB galaxies (jRieke et al.lll980h . experiencing a major SF 
formation episode in its nuclear regions, with strong super- wind and SN activity 



rormation episode m its nuclear regions, witn strong super-w ind and alN activity 
(Iwills et al.ll999h . Despite the many multi- wavelength studies ( Golla. Allen fc Kronberg 
19961 and references therein), only few measurements of nebular abundances, mainly 
derived from mid/far IR forbidden lines, have been published so far. From the anal- 
ysis of the [S ill] 35/i lin e in th e central 300 pc in radius at the distance of M82 
(~ 3.6 Mpc). lLorc et al. d200lh suggest an indicative [Si/H] ~ —0.4 dex relative to 
t he Solar value of Grevesse fc Saiival ( 1998h . 

iForster Schreiber et al.l (|200l| ). by analyzing the nuclear spectra of M82 taken 
with the Short Wavelength Spectrometer (SWS) onboard the Infrared Space Ob- 
servatory (ISO) within an aperture corresponding to 100 x 200 pc on the sky, give 
cold gas-phase abundances o f [Ne/Hl= +0.1, [Ar/Hl= +0.3 and [S/H]= —0.7 dex 
relative to the Solar values of Grevesse fc Sauval ( 1998| ). with an overall uncertainty 
between 0.2 and 0.3 dex. 

By taking the line ratios listed b y lAUoin. Co llin-Souffri n. fc Jolv which 
refer to the central ~ 100 pc of M82 (|Peimbert fc Spinradlll97ofr and the empirical 



calibrations of iMcGaughl (|1991[ ) which give N and Fe abundances as a function of 
[O/H] (see e.g. their Figure 11, equations 9 and 10a), one can also obtain rough 
estimates of [0/H]= +0.1, [N/H]= 0.0 and [Fe/H]= -0.3 dex relative to the Solar 
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values of Grevesse fc Sauval ( 1998f ). with an overall uncertainty of ±0.3 dex. 

The first atte mpt to measure ho t gas-phase abundances in M82 has been m ade 
quite recently by Ptak et aD ( 19971 ) and slightly revised by Umeda et al. (2002), by 
analyzing ASCA X-ray (0.5-10.0 keV) spectra. Fe and O abundances as low as 1/20 
Solar and significantly higher (by a factor between 3 and 1 0) Si, S, Mg, Ca ones 



have been obtai ned, also consistent with BeppoSAX results ( Gappi et al. 1999) 



More recently, Read fc Stevens! f)2002h by analyzing the spectra obtained with the 
Reflection Grating Spectrometer (RGS) onboard XMM-Newton found near-Solar 
Fe and O and super-Solar (between a factor of 2 and 5) Mg, Si, Ne, N abun- 
danc es. While the abu ndance ratios are in reasonable agreement within ±0.3 dex 



with iPtak et al.l ((1597), there is one order magnitude difference between the zero- 



point of the two calibrations. 



5.2 IR spectra and stellar abundances 



Near-IR spectroscopy is a fundamental tool to obtain accurate abundances of key 
elements like Fe, C, O and other metals (e.g. Si, Mg, Ca, Na, Al, Ti) in cool stars 
(X 



, g < 5000 K). Several atomic and molecular lines are stro ng, not affected by severe 



blending, hence als o measurable at low-medium resolution emmann Halllll 984 
lOriglia et all Il993t IWallace fc Hinklei Il997t I.Tovce et al.1 Il998t iFrogel et al.1 l200lh . 

making them powerful abundance tracers not only in stars but also in more dis- 
tant stellar clusters and galaxies and in a wide range of metallicities and ages 
( Origlia et al. 1997 : Oliva fc Origlia 1998^ 1 . Abundance analysis from integrated 



spectra of galaxies requires full spectral synthesis techniques to properly account for 
line blending and population synthesis to define the dominant contribution to the 
stellar luminosity. 

In SB galaxies the stellar IR continuum usually dominates over the possible gas 
and dust emission. This represents a major, conceptual simplification in population 
and spectral synthesis techniques, making possible and easier the interpretation of 
integrated spectra from distant stellar clusters and galaxies. 

Since the near IR c ontinuum of SB galaxies is dominated by luminous RSGs 
( Origlia fc Oliv^ll2000h . ;heir integrated spectrum can be modeled with an equiv- 
alent, average star, whose stellar parameters (temperature T e g, gravity Log g and 
microturbulence velocity £) mainly de pend on the stellar age and metallicity. Both 
observations and evolutionary models (|Kellerlll999t lOriglia et al.lll999t IOrigliall2003L 
and references therein) suggest that RSGs of ages between ~6 and 100 Myr and 
metallicities between 1/10 and Solar are characterized by low gravities (Log g < 1.0), 
low temperatures (< 4000 K) and relatively high microturbulence velocity (£ > 3 
km/s). A variation in the adopted stellar parameters for the average RSG popula- 
tion by AT cff ± 200 K, ALog g = ±0.5 and A£ =f 0.5 km s" 1 implies a < ±0.1 dex 
change in the abundances estimated from atomic lines, and < ±0.2 dex in those 
estimated from the molecular lines, which are more sensitive to stellar parameters. 

The 1.0-1.8// longslit spectra (see Fig. 15. lj) at R=2,500 of the nuclear region of 
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Figure 5.1: Near-IR spectra of the nuclear region of M82. Histograms: observed 
spectra; solid lines: synthetic stellar best-fit solution and two other models with ±0.3 
dex abundance variations. A few major stellar and nebular lines are also marked. 



M82 were obtained with the Near IR Camera Spectrometer (NICS) mounted at 
the Italian Telescopio Nazionale Galileo (TNG) on December 2002. The spectra 
were sky-subtracted, flat-fielded and corrected for atmospheric absorption using an 
O-star spectrum as reference. They were wavelength calibrated by using a Ne-Ar 
lamp and the mono dimensional spectra were extracted by summing over the central 
0.5" x 3.0", corresponding to an aperture projected on the source of 9 x 52 pc at the 
distance of M82. The spectra have been normalized to the continuum, which was 
determined applying a low-pass smoothing filter to each spectrum. Total integration 
times of 72 and 32 min in the J and H bands, respectively, were used, providing a 
final signal to noise ratio >40. 

By measuring the absorption line broadening a stellar velocity disp ersion a ~ 

'. 05=1=2 km/s has been derived, in perfect agreement with the estimate by Gafnev. Lester &: Telescol 
1993b and in the typical range of values measured in other massive SB galaxies 



Oliva et alJ ll999). 
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Figure 5.2: Left panel: probability of a random realization of our best-fit full 
spectral synthesis solution with varying the elemental abundances A[X/H] of ±0.2 
and 0.3 dex with respect to the best-fit (see Sect. 15. 2j) . Right panel: average difference 
between the model and the observed equivalent widths of a few selected lines (see 
Sect. O and Fig. EH- 



A grid of synthetic spectra of red supergiant stars for different inpu t atmospheric 



parameters and abundances have been com puted, using an upd ated (lOriglia et al 



120021 l2003h version of the code described in Origlia et al. ( 1993| ). Briefly, the code 
uses the LTE approximation and is bas e d on the molecular blanketed model at- 
mospheres of I Johnson. Bernat fc KnimJ (|l980h . It includes several thousands of 
near IR atomic lines and molecular roto- vibrational transitions due to CO, OH 
and CN. Three main compilations of atomic oscil lator strengths are u s ed, na mely 
the Kurucz's database 1 , a nd those published by iBiemont fc Orevessel l|l973h and 
Melendez fc Barbuvl (|l999h . 



The code provides full spectral synthesis over the 1-2.5 jum range. Given the high 
degree of line blending, the abundance estimates are mainly obtained by best-fitting 
the full observed spectrum and by measuring the equivalent widths of a few selected 
features (cf. Fig. 15. lj) . dominated by a specific chemical element, as a further cross- 
check. The equivalent widths have been measured by performing a Gaussian fit with 
a equal to the measured stellar velocity dispersion, typical values ranging between 
0.5 and 3 A, with a conservative error of ±200 mA to also account for a ±2% 
uncertainty in the continuum positioning. 

By fitting the full observed IR spectrum and by measuring the equivalent widths 
of selected lines, we obtained the following best-fit stellar parameters and abundance 
patterns for M82: T eff = 4000, Log g = 0.5, f = 3, [Fe/H]= -0.34; [0/Fe]= ±0.34, 
[<Si, Mg, Ca>/Fe]= ±0.38; [Al/Fe]= ±0.5; [C/Fe]= -0.26; 12 C/ 13 C<10. Table I5~TI 



1 ch. http : / / cf a-www. harward. edu/ amdata/ampdata/kurucz23/sekur .html 
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Table 5.1: Element abundances in Solar uni ts as derived from our ana lysis of the 
IR and X-ray spectra. The solar reference from Grevesse fc Sauval ( 1998| ) was used. 
Stellar abundances from TNG /NICS near IR absorption spectra of red supergiants 
(see Fig. 15.11 and Sect. 15.2)1 . Hot gas-phase abundances from XMM/RGS and pn 
spectra (see Figs. 15 A\ 15.61 and Sect. 15. 3j) . 



lists the derived abundances and thei r associated random error s at 90% confidence. 
Reference Solar abundances are from Grevesse fc Sauval ( 1998f ). 

Synthetic spectra with lower element abundances are systematically shallower than 
the best-fit solution, while the opposite occurs when higher abundances are adopted. 
In order to check the statistical significance of our best-fit solution, as a function 
of merit we adopt the difference between the model and the observed spectrum 
(hereafter 5). In order to quantify systematic discrepancies, this parameter is more 
powerful than the cl assical y 2 test, whi ch is instead equally sensitive to random and 
systematic scatters ( Qriglia et al. 20031 ). 

Since 5 is expected to follow a Gaussian distribution, we compute 5 and the 
corresponding standard deviation (a) for the best-fit solution and 6 test models with 
abundance variations A[X/H] = ±0.1, 0.2 and 0.3 dex with respect to the best- 
fit. We then extract 10000 random subsamples from each test model (assuming a 
Gaussian distribution) and we compute the probability P that a random realization 
of the data-points around a test model display a 5 that is compatible with the best- 
fit model. P ~ 1 indicates that the model is a good representation of the observed 



72 



5.3 X-ray spectra and hot gas-phase abundances 



spectrum. 

The left panel of Fig. 15.21 shows the average results for the observed J and H band 
spectra of M82. It can be easily appreciated that the best-fit solution provides in all 
cases a clear maximum in P (>99%) with respect to the test models. More relevant, 
test models with an abundance variation A[X/H] > ±0.2 dex lie at ~ 1.5a from the 
best-fit solution, while test models with A[X/H] > ±0.3 dex lie at o > 2 from the 
best-fit solution. 

The analysis of the line equivalent widths provide fully consistent results. The 
right panel of Fig. 15.21 shows the average difference between the model and the ob- 
served equivalent width measurements. Models with ±0.2 dex abundance variations 
from the best-fit solution are still acceptable at a ~ 1.5a significance level, while 
those with ±0.3 dex variations are only marginally acceptable at a 2 — 3<r level. 

Models with stellar parameters varying by AT e g ± 200 K, ALog g = ±0.5 and 
A£ =1=0.5 km s -1 and abundances varying accordingly by 0.1-0.2 dex, in order to still 
reproduce the deepness of the observed features, are also less statistical significant 
(on average only at > 2cr level) with respect to the best-fit solution. Hence, as 
a conservative estimate of the systematic error in the derived best-fit abundances, 
due to the residual uncertainty in the adopted stellar parameters, one can assume a 
value of < ±0.1 dex. 

By taking into account the overall uncertainty in the definition of the average 
population and the statistical significance of our spectral synthesis procedure, we can 
safely conclude that the stellar abundances can be constrained well within ±0.3 dex 
and their abundance ratios down to ~ 0.2 dex, since some (if not all) of the stellar 
parameter degeneracy is removed. 



5.3 X-ray spectra and hot gas-phase abundances 

The determination of the abundances of the hot X-ray emitting gas in SB galax- 
ies has traditionally suffered from large uncertainties. Indeed, the low angular and 
spectral resolution of the various X-ray telescopes in the pre-XMM-TVe wton/ Chandra 
era did not allow to disentangle between point sources a nd hot gas emission, mak- 
ing abundance determinations severely model-dependent (|Dahlem et al.ll2000h . The 



problem lies in the fact that the X-ray spectrum of SB galaxies contains two major 
components: the emission from hot diffuse gas and the integrated contribution of 
point sources. The first is described by an optically thin thermal spectrum plus emis- 
sion lines, the latter by a power-law. If the angular resolution is not good enough, 
it is not possible to reliably subtract the point sources from the total spectrum, so 
that the equivalent widths of the emission lines (and thus the element abundances) 
are not unambiguously defined. 

On t he other hand, high angular resolution alone is not enough: r ecent Chandra 
studies ( Strickland et al.ll2002i : Martin. Kobulnickv &: HeckmanlEopj) demonstrated 



that the relatively low spectral resolution of the ACIS detector makes individual 
element abundance analysis still problematic. The high spectral resolution of the 
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channel energy (keV) 

Figure 5.3: Total spectrum of the most luminous point sources in M82, obtained 
with Chandra. 



Chandra gratings rapidly degrades for sources which are much more extended than 
the instrumental Point Spread Function (PSF) and thus these gratings are almost 
useless for the study of nearby SB galaxies. However, by coupling the high angular 
resolution (0.5" PSF) of Chandra with the high spectral resolution (A/AA ~ 400) 
of the XMM/RGS, which is not sensitive to the source extension, it is possible to 
overcome the above described difficulties. 

M82 was observed several times by Chandra; we consider the two longest ex- 
posures, of 33 ks and 15 ks respectively, both obtained on September 1999. A 
cumulative spectrum of the brightest point sources was extracted from the two 
Chandra ACIS-I observations; the best-fit model was an absorbed power-law with 
N H = 7.9 ±0.7 x 10 21 cm" 2 , T = 0.84 ±0.07 x 2 r = 0-4, which is typical of High Mass 
X-ray Binaries (see Persic fc Rephaeli 20021 and reference therein). The 0.5-10 keV 
flux was 4.2 x 10 -12 erg s _1 cm -2 , which represented ~ 25% of the total (point 
source plus diffuse) flux (but n otice that M82 is a variable source, IPtak fc Griffith! 
1999L iReohaeli fc Grubeill2002l) . 

M82 was observed by XMM-Newton for 30 ks in May 2001; in the XMM-Newton 
observation, after light-curve cleaning, ~ 20 ks of data were available for scientific 
analysis. We extracted the EPIC pn spectrum from a circular region centered on 
the starbursting core of M82 and with a 15" radius, which approximately matches 
the RGS PSF; the standard recipe for EPIC spectra extraction suggested in the 
XMM web pages was followed. The XMMSAS software version 5.3.0 was used. 
A background spectrum was taken from standard background files provided by the 
XMM-Newton SOC. RGS spectra were extracted with standard settings (i.e. includ- 
ing events within 90% confidence and rejecting as background events those outside 
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Figure 5.4: XMM-RGS spectrum of the nuclear region of M82. Upper panel: 
data points (lower spectrum - open circles: RGS1; filled triangles: RGS2; upper 
spectrum: pn) along with the best-fit model. Lower panel: residuals in units of o. 
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95% o f the PSF figure); this procedure is similar to the one used by Read fc StevensI 
( 20021 ) but with more recent calibrations. After visual inspection, we discarded the 
2nd-order spectra because of low signal/noise ratio. We included in the fit only the 
lst-order channels comprised in the 6-20 A interval, since outside the S/N ratio 
rapidly drops. A systematic error up to 10% in the 14-18 A channels was included 
(XMM calibration document XMM-SOC-CAL-TN-0030 Issue 2). 

The hard-band (2-7 keV) Chandra images show that the point sources emission 
largely exceeds that of the hot gas, while the hard part (3-7 keV) of the pn spectrum 
can be well described by a power-law whose best-fit slope resulted in the same 
value previously determined for the Chandra point sources. Thus, making use of 
information from both Chandra and XMM-Newton observations, we have been able 
to constrain the point source slope, absorption and normalization at the moment of 
the XMM-Newton observation. 

We jointly fitted the pn and RGS spectra with XSPEC (Fig. 15 .4j) with a two 
component model which is the sum of: i) an absorbed power-law (accounting for 
point sources) with the slope and absorption parameters fixed at the best fit values 
obtained with Chandra while normalization was left free to account for variability; 
ii) an absorbed, optically-thin thermal plasma emission, with variable line intensities 
and differential emission measure (DEM) distribution. The DEM distribution was 

modeled with a 6^ order Chebyshev polynomial (c6pvmkl model in XSPEC); with 
the best-fit parameters it resembles a bell-shaped curve with a peak at kT ~ 0.7 
keV and FWHM ~ 0.65 keV ( i.e. the plasma tempera tures range from ~0.35 keV 
to ~1 keV), in agreement with lR.ead fc StevensI ()2002l ). The thermal component is 
absorbed by cold gas with best-fit column density Nh ~ 3.8 x 10 21 cm -2 . 

The confidence contours of the Iron abundance versus normalization of the ther- 
mal component are reported in Fig. 15.51 As expected the two parameters are well 
correlated as the higher is the normalization of the continuum spectrum the lower 
is the intensity of the emission lines and in turn the element abundances. The 
derived abundances and their associated 90% confidence error are shown in Ta- 
ble 15.11 We notice that th ese abundances are about a factor 2 below those reported 
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R,ead fc Steve^l (|2002lk we attribute the discr epancy to the different modeling 



of the continuum emission fe.ead fc StevensI 120021 did not include pn data, neither 
considered the contribution from point sources). However, it is worth of mentioning 
that both sets of abundances indicate a low O abundance compared to the other 
a-elements. 

To further check the consistency and robustness of our estimates, we repeated 
the fit two times considering only the RGS data and using i) the above described 
thermal plus power-law model and ii) a thermal model alone. T his last check, in 
partic ular, most closely matches the analysis previously made by Read fc StevensI 
(2002) to the point that the discrepancies are within the errors and can be attributed 
to the differences in the used instrument calibration. The abundances derived from 
the three different models are plotted in Fig. 15.61 Those obtained combining pn and 
RGS have significantly smaller uncertainties. One can see that different models have 
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Figure 5.5: 68, 90 and 99% confidence contours of the Iron abundance vs in- 
tensity of the thermal component. For ease of reading the y-axis is shown both in 
linear and logarithmic scale. 
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Figure 5.6: Best-fit element abundances of the hot gas for different models. 
Filled circles: fit of both RGS and pn data, thermal plus power-law model. Open 
squares: RGS data alone, thermal plus power-law model. Filled triangles: RGS data 
alone, thermal model alone. For ease of reading the y-axis is shown both in linear 
and logarithmic scale. 
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Figure 5.7: Hot gas-phase Oxygen abundance vs absorption. For ease of reading 
the y-axis is shown both in linear and logarithmic scale. The correlation, if any, is 
rather weak, and the Oxygen abundance can be well constrained. 

mainly the effect to almost rigidly scale all the abundances, with minor impact on 
the abundance ratios. 

The low Oxygen abundance is somehow difficult to explain in the framework of 
a-element enhancement by type II SN explosions. Thus, to check its reliability 
we included in the fit also the channels corresponding to the 20-23 A interval, 
where the Ovil triplet is observed (although with very few counts). We find no 
significant difference neither in Oxygen abundance nor in DEM distribution (the 
O VIII A19/0 VII A22 ratio is extremely sensitive to the plasma temperature). Ne is 
also somewhat under- abundant. We checked for possible instrumental effects, such 
as variations in the RGS effective area in proximity of the Ne and O lines; however 
the effective area in this region is a rather smooth function of the wavelength, so 
that there should be no instrumental issues and these low abundances are likely to 
be real. 

Since M82 is observed edge-on, the X-ray emission of its nuclear region is subject 
to quite heavy absorption which could affect the O lines detected at low energies. 
However, as shown in Fig. 15.71 the contour plot of the Oxygen abundance towards 
the foreground absorption, suggest a rather weak correlation. In the next section 
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Figure 5.8: Stellar, hot and cold gas abundance ratios relative to Iron. Cold 
gas-phase abundances are from different data sets in the literature, hence they are 
not homogeneously determined. The dashed line indicates the solar values. For ease 
of reading the y-axis is shown both in linear and logarithmic scale. 

we discuss other physical possibilities to explain the O under-abundance issue. 



5.4 Discussion 

In M82 all the three components, namely hot and cold gas-phases and stars trace a 
very similar Iron abundance, the average value being [Fe/H]~ —0.35 dex. Indeed, 
since Iron is mainly produced by SNe la, it is expected to be released in the ISM 
only after ~ 1 Gyr from the local onset of SF. 

At variance, a— elements (O, Ne, Mg, Si, S, Ca, Ti) are predominantly released by 
SNe II with massive progenitors on much shorter timescales. Stars trace an average 
[<Si,Mg,Ca>/Fe]~ 0.4 dex, while the hot gas suggests an average [<Si,Mg,S>/Fe]~ 
0.5 dex. Such an overall [a/Fe] enhancement in M82 is fully consistent with a stan- 
dard chemical evolution scenario, where the ISM in the nuclear regions of massive 



SB g alaxies is mainly enriched by the products of SNe II explosions (see e.g. lArnett 



19951 for a review) occurring in recursive bursts of SF of relatively short duration. 
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Fig. 15.81 shows the various abundance ratios relative to Iron, as measured in stars, 
hot and cold gas-phases. Since the cold gas-phase abundances are from different data 
sets in the literature, hence not homogeneously determined, the inferred abundance 
ratios should be regarded with caution. 

When comparing the hot and cold gas-phase abundances, one finds Ne over- 



abund ant, Si and especially S significantly under-abundant in the cold gas. iForster Schreiber et al 



(2001) suggests S and to a lower level Si depletion onto interstellar dust grains. How- 
ever, other metals like Fe in par ticular, with even higher degree of incorporation into 
dust ( Savage fe Sembach 19961 ). should be severely depleted in the cold gas, which 
does not seem the case . A possible explanation is that the o ptical lines used to infer 



Fe and O abundances (jAUoin. Collin- Souffrin. fc .T olv 
lines used to infer Ne and S abundances ((Lord et al 



1979) and the the mid /far IR 



2001; Forster Schreiber et al. 



20011 ) trace different nebular sub-structures within the central few hundreds pc, 



somewhat chemically dishomogeneous or with different dust content. 

The [O/Fe] abundance ratio in M82 is even more puzzling to interpret. Before 
the onset of the current burst of SF, the ISM was enriched in O as well as in the 
other a-elements, as suggested by our stellar and cold gas-phase abundances. The 
O under- abundance measured in the hot gas cannot be modeled with a sta ndard 
nucleosynthesis from the present generation of SNe II. lUmeda et al. (|2002l) sug- 
gested explosive nucleosynthesis in core-collapse hypernovae but they predicted a 
large under- abundance of Ne and Mg as well, which is not observed. [O/Fe] under- 
abundance can be also explained with a major Fe enrichment by current SN-Ia explo- 
sions from previous generations of stars, but in this case all the [a/Fe] ratios should 
be low, again contrary to what observed. O might be also locked into dust grains. 
Indee d, SNe II are known to produce significant amounts of dust (|Todini fe Ferrara 



20011 ) and O, being a basic constituent of dust grains, is expected to suffer of severe 
dust depletion. However, ot her metals like Fe in part icular, but also Si and Mg 
should be depleted into dust ( Savage fe Sembach 1996| ). leaving almost unchanged 
or eventually further enhancing the [a/Fe] ratio. 
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6 Summary 



We have analyzed a small, but well denned sample of 17 star-forming galaxies, 
extracted from the HFS97 catalogue, for which there is a homogeneous information 
on op tical, FIR, radio an d X-ray bands (local sample). In agreement with previous 
work ( David et al.lll992l ) we find that the logarithms of the soft (0.5-2 keV) X-ray 
luminosities (corrected for Galactic absorption only) are linearly correlated with the 
logarithms of both radio (1.4 GHz) and FIR luminosities. We have extended our 
analysis to the harder X-ray band, essentially free from internal absorption which 
may affect the soft X-ray fluxes, and found that there is a tight linear correlation 
between the X-ray luminosities in the 2-10 keV interval with both the radio and the 
FIR luminosities, normally assumed as the indicators of the star formation rate. The 
addition of 6 galaxies (supplementary sample) homogeneous with, but not included 
in HFS97, does not modify these results. We conclude that the origin of the hard 
X-ray emission must be closely related to star formation and calibrate an X-ray SFR 
indicator. 

Candidate starburst galaxies have been selected in the HDF North, with redshifts 
up to z ~ 1.3, and their rest-frame X-ray luminosities are computed by extracting 
counts in redshifted bands from the Chandra observation of the HDFN. With this 
approach we have shown that the 2-10 keV/radio linear correlation holds up to 
z ~ 1.3, encompassing five orders of magnitude in luminosity, up to L2-10 ~ a f ew 
10 42 erg s -1 and a corresponding star formation rate ~ 1000 M© yr" 1 — but at 
such large luminosities some contribution from nuclear activity might be present. 
The fit to the 0.5-2.0 keV/radio data is also linear up to z ~ 1.3. 

As an additional investigation we have also analyzed a sa mple of LLAGN (LINERs 
and Seyfert's) included in HFS97 ( Terashima et al. 2002| ): while, as expected, the 
X-ray luminosities are generally in excess with respect to star-forming galaxies for 
the same FIR luminosity, the distribution of the objects in the X-ray vs. FIR lumi- 
nosity diagram is bounded from below from the region occupied by the star-forming 
galaxies, indicating that the X-ray emission of LLAGN falling in this border-line 
region could be mainly due to star formation processes, rather than being of nuclear 
origin. 

We used the radio/FIR/X-ray correlations to convert the local (z = 0) FIR and 
radio luminosity functi ons (LFs) into X-ray LFs, which were compared with the 
one recently derived in iNorman et al. ( 2004 ) for objects at intermediate redshifts 
(0.1 < z < 1.2). A quite good agreement was found for L x < 10 41 erg/s, while 
at larger luminosities there is almost an order of magnitude of difference. Tho 
possible explanations may be proposed: i) a fraction of AGN is contaminating the 
high-luminosity tail of the sample of galaxies used in deriving the X-ray LF; ii) the 
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infrared and radio surveys do not sample a sufficiently large volume to find rare, 
high luminosity objects. 

X-ray number counts for the normal galaxies were predicted both by integrating 
the radio- and FIR-based luminosity functions, and by converting the radio Log N- 
Log S for the sub-mJy population. The expected counts extend much below the 
limiting fluxes of the deepest X-ray surveys (about one order of magnitude in the 
soft band, one and a half for the hard band), a nd are within the limits s et by the 
fluctuations analysis in the Chandra deep fields (|Mivaii fc Griffithsll2002al lbh. They 



are also co nsistent wit h the p redictions based on the evolution of the cosmic SFR 



density bv iPtak et al.l ([20011 ). Since the results from fluctuations analysis in the 
soft X-rays suggest an excess of sources with respect to AGN synthesis models at 
fluxes below ~ 10~ 17 erg s _1 cm -2 , it may be possible that the sub-mJy galaxies 
represent the dominant population in the X-rays at very faint fluxes. Our analysis 
of the X-ray number counts also allows a determination of the fraction of normal 
galaxies among the X-ray sources detected in the deepest X-ray surveys. Namely, 
we estimate that 20% of the X-ray sources with fluxes larger than 5 ■ 10~ 17 erg s _1 
cm -2 in the 0.5-2.0 keV band are star forming galaxies. 

The contribution to the cosmic X-ray background in the 2-10 keV band is esti- 
mated by integration of the derived X-ray number counts. The contribution from 
galaxies detected in the deepest radio surveys is < 2%. This estimate may rise up 
to 11% by extrapolating the radio counts down to 1 /xJy, or ~ 10~ 18 erg s _1 cm -2 
in the X-ray band. However, since a fraction of the sub-mJy objects may not be 
star-forming galaxies, these figures for the time being should be regarded as upper 
limits. 

Making use of our calibration of the X-ray luminosity as a SFR indicator, it was 
also possible to rescale the X-ray luminosity function into a cosmic star formation 
rate density. This definitely assesses the possibility of using X-ray observations to 
obtain an absorption-free estimate of the cosmic star formation history. However, 
the deepest surveys available do not detect normal galaxies beyond z ~ 1.3; the 
most interesting period of the cosmic history (from z ~ 2 to z ~ 4, where it is 
still debated if the cosmic star formation rate density declined or stayed constant) 
will remain precluded to study until future X-ray telescopes, such as XEUS, will be 
launched. 

Our abundance analysis of the stellar and hot gas-phase components in the nu- 
clear region of M82 indicate an Iron abundance about half Solar and an overall 
a— enhancement by a factor between 2 and 3 with respect to Iron. These abun- 
dance patterns can be easily explained within a standard nucleosynthesis scenario 
where the ISM is mainly enriched by SNe II on relative short timescales and with a 
star formation process occurring in recursive bursts. 

Oxygen behaves in a strange fashion. It is over-abundant in stars and cold gas, 
similarly to the other a-elements, while is significantly under-abundant in the hot 
gas. Major calibrations and/or modeling problems seem unlikely. Hypernovae nu- 
cleosynthesis, dust depletion, SN-Ia enrichment can somehow explain an Oxygen 
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under-abundance but other metals should follow the same pattern, while they do 
not. Somewhat exotic threshold effects and/or depletion mechanisms preferentially 
affecting Oxygen could be at work in the nuclear region of M82, but presently this 
issue remains controversial. The situation will be probably clarified in the next fu- 
ture, since we were recently granted a five-times-deeper reobservation of M82 with 
XMM- Newton. 
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Plate I: Spitzer observation of M81. The 3.6/i emission is mainly due to red super- 
giant stars, while the emission at longer wavelenghts is progressively dominated by 
dust in H II regions. 
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Plate II: Spitzer observation of M81. The 24// and 70/i emission is mainly due to 
warm (~ 40K) dust in H II regions, while the 160/z emission traces a colder (~ 20K) 
dust component. 
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(a) X-ray emission from the Antennae, showing the low angular resolution of past 
X-ray telescopes. Right panel: optical image from the Digitized Sky Survey, with 
countours superimposed. Left panel: ROSAT PSPC image with optical contours 
superimposed. Middle panel: ASCA GIS image with optical contours 
superimposed. 




Plate III: (b) The Antennae galaxy observed with Chandra, from iFabbiano et al 



(200^). In the upper and lower left panels the colour coding represents the energy of 
the incoming photons (red: 0.5-1 keV; green: 1-2 keV; blue: 2-7 keV). In the lower 
right panel, the color coding is referred to the different chemical elements present in 
the spatially resolved spectrum (red: Fe; green: Mg; blue: Si). 
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Plate IV: Chandra image of M82. The colour coding represents the energy of the 
incoming photons (red: 0.5-1 keV; green: 1-2 keV; blue: 2-7 keV). The field size of 
the upper panel is ~ 15' x 15'. 
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Plate V: Chandra (upper panel) and XMM-Newton (lower panes) images of M82. 
The (smoothed) image in the upper panel shown the inner region of M82. The 
lower panel has a field diameter of ~ 20' and shows the different angular resolution 
of XMM-Newton (note that the image is shown with a different position angle). The 
colour coding represents the energy of the incoming photons (red: 0.5-1 keV; green: 
1-2 keV; blue: 2-7 keV). 
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(a) XMM- Newton image of NGC253. The central starburst is clearly visible (the 
large white spot at the centre of the image). Also the galactic disc is visible (red 
diffuse emission). Colour coding as in the previous plates. 
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Plate VI: (b) Chandra images of NGC253, from IStrickland et all (|200nh . showing 
the different morphology of the soft and hard X-ray emission. The galaxy outflow 
is clearly visible in the left panel, while the hard X-ray emission is dominated by 
sources concentrated in the galaxy centre. 
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Plate VII: True color, smoothed Chandra image of the centre of NGC 3256 (the 
actual colour is the sum of red, green and blue, the intensity of each one representing 
the flux in the 0.3-1.0 keV, 1.0-2.0 keV and 2.0-8 keV respectively); north is up, 
east is left; the distance between the two cores is ~ 5". The two cores are the 
brightest sources both at 1.4 GHz and in the 2-10 keV band, but since the southern 
one lies behind a dust lane, its soft X-ray emission is completely absorbed. While 
the northern core follows both the radio/soft X-ray and radio/hard X-ray relations 
(Eq. I2.fij2.10ft . the southern one only follows the radio/hard X-ray relation. 
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